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Abstract. Strichartz-type estimates for one-dimensional surface water-waves 
under surface tension are studied, based on the formulation of the problem as 
a nonlinear dispersive equation. We establish a family of dispersion estimates 
on time scales depending on the size of the frequencies. We infer that a solu- 
tion u of the dispersive equation we introduce satisfies local-in-timc Strichartz 
estimates with loss in derivative: 

2 1 _ 1 

II"IIlp([0,T])W=-1/p.<!(M) ^ „ ~ o' 

where C depends on T and on the norms of the initial data in X H''~^^^. 
The proof uses the frequency analysis and scmiclassical Strichartz estimates 
for the linealized water-wave operator. 
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1. Introduction 

The problem of surface water waves, in its simplest form, concerns the two- 
dimensional dynamics of an incompressible inviscid liquid of infinite depth and the 
wave motion on its one-dimcnsional surface layer, under the influence of gravity and 
surface tension. The moving surface is given as a nonsclf-intersecting parametrized 
curve. The liquid occupies the domain below the curve, where the liquid motion 
is described by the Eulcr equations under gravity. The flow beneath the moving 
surface is required to be irrotational. The kinematic and dynamic boundary condi- 
tions hold at the moving surface, stating respectively that the normal component 
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of velocity is continuous along the moving surface and that the jump in pressure 
across the moving surface is proportional to its mean curvature. The flow is as- 
sumed to be almost at rest at great depths, and the moving surface is taken to be 
asymptotically flat. 

Provided with the initial surface profile and the initial state of fluid current, the 
water-wave problem naturally poses as an initial value problem. Early mathemat- 
ical results for local well-posedness date back to [17,23] and include [12,22,33,34]. 
Following the works by Sijue Wu [30, 31] there has been considerable progress 
in the study of local well-posedness for the water-wave problem as well as for a 
class of the Euler equations with free boundary. We refer the interested reader 
to [3,4,10,11,20,21,25], and references therein. Recently, results for long-time 
existence [24, 32] appeared for gravity water waves of infinite depth. 

Nonlinearity characteristic to the boundary conditions at the moving surface 
significantly restricts the range of analytical tools available for the existence theory 
for the water-wave problem. As a matter of fact, all results listed in the previous 
paragraph on local well-posedness hinge upon obtaining high energy expressions 
and establishing their bounds, namely the energy method. While construction of 
such energy expressions is nontrivial and design of an iteration scheme is involved, 
nevertheless, results from the energy method do not provide any further information 
about solutions, other than that they remain as smooth as their initial states. Better 
understanding of the dynamics of surface water waves can be made with the help 
of a priori estimates other than energy estimates. 

On the other hand, the dispersion relation (see Remark I2.6P 



of surface water waves provides a guiding principle of their linear dynamics. Here, 
c{k) is the speed of the simple harmonic oscillation with the wave number k; S ^ 
is the coefficient of surface tension and g ^ is the gravitational constant of 
acceleration. Under the influence of surface tension, i.e. S > 0, the fact that the 
phase velocity c{k) is asymptotically proportional to the square root of fc as fc — > oo 
indicates a certain "regularizing" effect by the process of broadening out the surface 
profile. In the gravity-wave setting, i.e. 5 = and g > 0, in contrast, p.ip does 
not induce such a regularizing effeci|j. 

Dispersive properties have paramount importance in the study of nonlinear 
Schrodinger equations, the Korteweg-de Vries equation, nonlinear wave equations, 
and other nonlinear dispersive equations. In the recent work of Alazard, Burq and 
Zuilly [1], local smoothing effects are obtained for water waves under surface ten- 
sion (see also Appendix [XJ. Such a smoothing effect is a direct consequence of the 
dispersive property of surface water waves, and it contrasts markedly with what 
energy estimates alone can tell. The present purpose is to investigate the dispersive 
property for the water-wave problem with one-dimensional surface under surface 
tension by establishing estimates of the solution imder the mixed Sobolev norms, 
commonly referred to as Strichartz estimates. 



Gravity waves may still be thought of as "dispersive" in the sense that wave components 
with different frequencies propagate at different speeds; see [32], 



(1.1) 
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1.1. The main results. The present treatment of the dispersive property for the 
water-wave problem under surface tension is based on the formulation of the prob- 
lem as a second-order in time nonlinear dispersive equation 

(1.2) d'^u - ^Hdlu + gHdaU = ~2udtdaU - u^dlu + R(u, dtu), 

coupled with a transport-type equation (|3.8p . We shall derive it in Section [2] and 
Section [3] Here, u is related to the tangential velocity at the moving surface and 
it serves as the unknown; t £ M_|_ is the temporal variable and a € R is the 
(renormalized) arclength parametrization of the curve, which serves as the spatial 
variable. The Hilbert transform, denoted by H, may be defined via the Fourier 
transform as Hf{£^) = — isgn(^)/(^). The remainder R is of lower order compared 
to 2udtdaU and u'^d'^u in the sense that 

\\R{u,dtu)\\Hs C{\\u\\H.+i,\\dtu\\H^) 

for s ^ 1. Here and elsewhere, H'^ means the L^-Sobolev space of order s in the 
variable a g M. 

Our main results concern Strichartz estimates for the water-wave problem under 
surface tension with loss in derivative. In the course of the proof, its local wcU- 
posedncss is proved. 

Theorem 1.1. Let S > and g ^ be held fixed. For s > 2 + 1/2 the initial value 
problem of p.2p prescribed with the initial conditions 

u(0,a) — uo{a) and dtu(0, a) — ui{a), 

where {uo,ui) G i7'*(IR) xi7''~'^/^(M) is locally well-posed on a time interval t G [0,T] 
for some T > 0, and the solution u satisfies {u{t) , dtu(t)) G C([0, T]; iJ*(]R) x 

H''-3/2(K)). 

Moreover, if s is sufficiently large, the solution u satisfies the inequality 

(1.3) (y (y \D'-^/'Pu{t,a)\Ua) dtj ^ C, 

where (p, q) satisfies the admissibility condition 

, , 2 11 

1.4 - + -=o, g<oo, 

p q 2 

and C > depends on s,q,p,T and \\ui)\\u^(K)i \Wi\\h'-^/^(k.)- Here and in sequel, 

Da = -ida- 

If the solution is localized to dyadic frequency bands and semiclassical time 
scales, the estimate is better. 

Theorem 1.2. Let %jj^{Da) be a Fourier multiplier supported in frequencies 2^~^ ^ 
1^1 ^ 2-'"'"^. Under the hypothesis of Theorem \1.1\ with s sufficiently large, the 
frequency-localized solution 'ip^(Da)u satisfies 

(1.5) (y (y \D'-^/^^i,^{Da)u{t,a)\'ida) dtj ^ C, 

where {p,q) satisfies (|1.4p with g ^ oo and C > depends on s,q,p,T and 

II'"o||h^, ||ui||^.-3/2. 
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Notations. Recorded here are the notations and conventions used in the sequel. 

Let ^ k,l ^ oo and 1 ^ p, g ^ oo. By W^''*(]R) we mean the L"^ Sobolev space 
on a G M of order k, and by W;'P{[0,T]) we mean the Sobolev space on the 
interval t S [0,T] of order I. By Hl{[0,T]) the Sobolev space on the interval 
t G [0, T] of order / . We will also use the Sobolev spaces of negative order, Hl^{R) 
with fc < 0. For ^ p, g ^ oo we recall the definitions for the mixed Sobolev spaces 
i«(R)Lf([0,r]) and ([0, T])L« (M) by the norms of these spaces 

II/|IlJ(R)LP([0,T]) 



\\f\\miO,T])Ll{R) 

We write L^L^, for X« (M)if ([0, T]) and L^i« for if ([0, r])L9 (R) when there 
is no ambiguity. We use the analogous convention for W^-'^W^^'^, W^^W^^''^, and 

I. 2. Perspectives. The derivative loss of 1/p in Theorem 11.11 is likely not sharp, 
as the following heuristic arguments indicate. 

For any dispersive equation in one spatial dimension, the W^'^ decay 
rate is t~^^^ (with loss of s derivatives depending on the equation). If we linearize 
about the zero solution (see (|1.14p below), we see the solution satisfies Strichartz 
estimates with the admissibility condition (|1.4p and a l/2p derivative loss (see, for 
example, [8]), which is an improvement of l/2p derivatives compared to Theorem 

II. 11 Moreover, this equation satisfies the scaling symmetrjH 

u{t,a) ^ X^^^u{X^^^t, Xa) 

for any dilation factor A > 0, and we readily verify that Strichartz estimates with 
the admissibility condition (|1.4p and l/2p derivative loss is invariant with respect 
to this scaling. We thus expect Strichartz estimates with admissibility condition 
(jl.4p and l/2p derivative loss to be optimal. That is, Theorem ll.il represents twice 
the loss in derivative of the optimal estimate. 

However, this optimal estimate cannot be obtained by interpolation with known 
estimates, even in weighted form. Indeed, to compare to the local smoothing esti- 
mate ( [1] or Appendix [X)) . if we use Sobolev embeddings, we have 

and if we use that Dt is comparable to Da^'^ (at least for a solution linearized about 
0), in turn, we have 



^In the absence of the effect of gravity, g = 0, the nonhnear equation 1 11.211 also enjoys this 
scaling symmetry. This follows from the scaling symmetry of the Euler equations and the dynamic 
boundary condition that the jump of pressure across the moving surface is proportional to the 
mean curvature of the surface. 



<?/p 



1/9 



\f{t,a)\Pdt\ da 



|/(t,a)|«da) dt] 
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Suggested by scaling 

Sobolev embedding plus local smoothing 

Holder plus energy 




1/2 



1 



l/p 



Figure 1. Fixed time scale, l/2p derivative loss. The {p,q) re- 
lation suggested by scaling, from Sobolev embeddings plus local 
smoothing effect, and from Holder's inequality in time with Sobolev 
embeddings and energy conservation. 



By the local smoothing effect gain of 1 /4 derivative, we then bound || (a) Da ^^^''uj 
p > 1/2, in terms of the initial data in iJ''(R) x 7J*"^/^(R), provided that 

2 1 

- + - = 1. 

p q 

This is weaker than the optimal estimate. On the other hand, if we use Holder's 
inequality plus energy conservation, we get a loss of l/2p derivatives provided 

1 1 _ 1 
2^ ^ ^ ^ 2 

(see Figured]). 

To make a direct comparison of the estimate of Theorem 11.11 with the optimal 
condition is not as clear, since we must use Sobolev embeddings somewhere. If 
we do use an additional Sobolev embedding in the discussion above to make a 
comparison of l/p derivative loss, the optimal admissibility condition becomes 

5 1 _ 1 
2^ ^ 9 ^ 2' 

while that from smoothing is (|1.6p with the right hand side replaced by 1, and that 
for energy estimates is 



(1.6) 



1 
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1 



Theorem 11.11 

Suggested by scaling and Sobolev 

Sobolev embedding plus local smoothing 

Holder plus energy 



1/2 





Figure 2. Fixed time scale, 1/p derivative loss. The {p,q) rela- 
tion given in Theorem ll.il that suggested by scaling, from Sobolev 
embeddings plus local smoothing effect, and from Holder's inequal- 
ity in time with Sobolev embeddings and energy conservation. 



(see Figure [2]). Again we see that the estimate of Theorem 11.11 cannot be obtained 
by interpolation between known estimates. 

On the semiclassical time scale ^ < ^ 2^-'/^T, our Strichartz estimate (|1.5p has 
a smaller loss in derivative, and the optimal scaling condition is the same as (jl.4p . 
Since the local smoothing cannot be improved on the semiclassical time scale, our 
estimate ()1.5p represents a larger gain over what Sobolev embeddings plus local 
smoothing could tell us on the semiclassical time scale (see Figure [3]) . 

1.3. Idea of the proofs. While (|1.2p is dispersive, its nonlinearity is severe, and 
as such in the study of its dispersive properties one must take its nonlinear effect 
into account. To better understand the strength of nonlinearity versus the weak- 
ness of dispersion we examine the local smoothing effect for (|1.2p . An application 
of Parseval's formula, together with a change of variables, shows that ( [19] for 
instance) the solution of the linear homogeneous equation 



gains 1/4 derivative of smoothness over the initial data. An application of a TT* ar- 
gument then shows that the solution of the corresponding inhomogeneous equation 
gains 2 derivatives of smoothness over the inhomogcneity. But, this local smoothing 
effect is not enough to control nonlinear terms in (|1.2p containing more than two 
spatial derivatives, e.g. 2udtdaU. 



(1.7) 




u = 



S>0 
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1/9 



Theorem 11.21 (and scaling) 

Sobolev embedding plus local smoothing 

Holder plus energy 




Figure 3. Semiclassical time scale, l/2p derivative loss. The 
(p, q) relation given in Theorem 11.21 fagrees with that suggested by 
scaling), Sobolev embeddings plus local smoothing, and Holder's 
inequality in time plus Sobolev embeddings plus energy conserva- 
tion. 



To overcome this setback and to obtain Strichartz estimates for the solution of 
nonlinear equation p.2p . we write it as 

(1.8) d'^u - -Hdlu + gHdaU + 2udtdaU + u^dlu = R[u, dtu). 

That is, we view 2udtdaU and u^d^u as "linear" components of the equation, but 
with variable coefficients which happen to depend on the solution itself. In other 
words, we reduce the size of nonlinearity at the expense of making its linear part 
more complicated. We then make a serious effort to establish Strichartz estimates 
for the linear operator 

(1.9) df - ^Hdl + gHd^ + 2V{t, a)dtda. + V\t, a)dl, 

for a class of functions for the variable coefficient V{t, a). 

The operator p.9p may be thought of the operator — Hd^ perturbed by 
variable-coefficient but lower-order terms 2V{t^a)dadt + V^{t^a)d^. While the 
added terms are of lower order they are not constant, and they bring a great deal 
of difficulty in the analysis of the paper, which is the heart of the matter. 

In [1] and in Appendix lA} in order to establish the local smoothing effect for the 
nonlinear equation (|1.2p . similar approaches are employed. 
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1.3.1. Construction of the parametrix. Our approach to cstabhshmg microlocal 
Strichartz estimates for ()1.9p is based on the construction of its approximate solu- 
tion. 

When V{t, a) = 0, the solution of the homogeneous equation (|1.7p is given by 
the formula 



where uq and ui describe the initial data. Here, for the sake of exposition, we have 
assumed S/2 ~ 1 and g = 0. Motivated by this, we make an oscillatory integral 
ansatz 



to solve the problem associated to (|1.9p . The phase functions ip^ is chosen to satisfy 
(y5^(0, a,^) = a^, and as such the recovery of the initial conditions entails solving 
for a system of elliptic pseudodifferential equations. 

Applying the linear operator (|1.9p to our ansatz, we consider the worst terms, 
produced when first-order derivatives fall on the phase functions. They make a 
first-order nonlinear equation (|4.12p for ip^, commonly referred to as the eikonal 
or Hamilton- J acohi equation. The usual approach to solving the Hamilton- Jacobi 
equation is through the technique of generating functions for the associated Hamil- 
tonian. The equation (|4.12p is, however, neither homogeneous nor polyhomoge- 
ncous (in tp^ and ip^)., and as such solutions are found on a time scale comparable 
to See Lemma l475l for details. We thus construct phase functions for each 

dyadic frequency band |^| ^ 2^ on a frequency-dependent time scale t ^ 2~^/^. The 
construction of the leading- order parametrix w is detailed in Section |4l 

1.3.2. Semiclassical Strichartz estimates. We explain our strategy to establish Strichartz 
estimates for the linearized water-wave operator (|1.9p under surface tension. 

Let us first discuss basic ideas for Strichartz estimates for the one-dimensional 
free Schrodinger equation 



since we will use similar ideas. Prescribed with the initial condition w(0, a) = UQ{a), 
the solution of (jl.lip can be written via the Fourier transform as 



(1.10) 





(1.11) 



idtu + dlu = 0, 




We write this as a convolution with an integral kernel as 





The phase function <^(^; t, a, 0) = ^(a — /3) + t^ 



^ has a critical point at 



^^(^(^e) = a - /3 + 2i^e = 0, or ^, = {P-a)/2t. 
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Since d'^(p{^c) ~ 2t, moreover, the phase is nondegenerate for t > 0. Then, by the 
standard method of stationary phas^l, at least for uq locahzed in frequency, we 
obtain K = Ki + (smoothing) with 

(1.12) \K,{t,a,p)\i:cr^/^ 

where C > is independent of i, a and /3. 

Next, we recall an abstract result which follows from the work of Ginibrc-Vclo [15, 
16] and recorded in the paper of Keel and Tao [18] , stating that a dispersion estimate 
leads to Strichartz estimates under the L^-norm for a range of (p, q) depending 
on the strength of the dispersion (the power of t in the dispersion estimate). 

Theorem 1.3. Let (X,dx) be a measure space, let H be a Hilbert space, and let 
U{t) : H L'^{X) be a linear operator satisfying 

(i) \\U{t)f\\Li^C4f\\H, and 

(ii) \\U{t')U*{t)gU^^C,\t-tT\\gUi 
for some cr > 0. Then for every pair (p, q) satisfying 

1 a a 

p^l^r 

the estimate 

\\U{t)f\\L'^Ll^C2\\f\\H 

holds true, where C2 > depends only on Ci, a, p and q. 

The semiclassical dispersion estimate wc prove in this paper depends also on the 
semiclassical parameter . A rescaling in time and application of Theorem 11.31 
gives the following semiclassical Strichartz estimate theorem (see, for example, [14, 
Theorem B.IO]). 

Theorem 1.4 (Semiclassical Strichartz estimates). Let {X, dx) be a measure space, 
let Hq > fixed, let H be a Hilbert space, and let U(t) : H L'^{X) be a linear 
operator satisfying 

(i) \\U{t)f\\Li^C4f\\H, and 

(ii) \\U{t')U*{t)g\\L^^C2h-^^\t-t'\-''\\g\\Li 

for some a > and all < h ^ ha. Then for every pair (p, q) satisfying 

1 (T (T 

p + g = 2' 

the estimate 

\\Uit)f\\^Li^C3h-i^\\f\\H 

holds true, where C3 > depends only on Ci, C2, o, [i, p and q. 

In light of the above theorem, (|1.12p gives that a solution of (|l.lip satisfies the 
estimate 

hh-Ll ^ C\\uo\\l- 

where {p, q) satisfies 

2 11 

p + g " 2- 

Estimate 111.121 1 is usually derived from the explicit formula for the kernel K, but here we want 
to stress a method that can be generalized for variable coefficient dispersive differential operators. 
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Furthermore, a scahng argument assures that the estimate is sharp. 

Returning to our setting, we consider the hnear problem 

(dfU - HdlU + 2V{t, a)da,dtU + V^{t, a)dlU = R{t, a) 
([/(0,a) = C/o(a) and dtU{Q,a) = Ui{a), 



(1.13) 



on the time scale [0, 2~-'/^T]0 where U,R,Uo and Ui are localized to the dyadic 
frequency band 2^~^ ^ |^| ^ 2''+^. Here, for simplicity we take S/2 = 1 and 5 = 0. 
We write the oscillatory integrals 



(1.14) j j e±^*l«l'''e^«("-'^'(f^o(/3) T i|er^/';7i(/3))dM 

associated to the solution of the zero-coefficient equation, V{t,a) = 0. 

Considering the corresponding phase for ^ large and positivqj, t > 0, and with 
the + sign, let ip{£,; t, a, (3) — ^(a — /?) + i^^^^. Its critical point is at 



(1.15) d^v{i,)^-tej^ + a-l3 = Q, or g 



Since 



the critical point is nondegcnerate for t > and we are ready to use the method 
of stationary phase. However, plugging these results into the stationary phase 
argument docs not yield a bound on the kernel uniform in a or (3. 

To remedy this, we use propagation of singularities to estimate {a — f3)/t in 
terms of derivatives in (3 to obtain a dispersion rate of t~^/^ with loss in derivative. 
But, our parametrix is for the operator (|1.9p with variable coefficients, which is 
considerably more complicated than (|1.14p . Moreover, the parametrix exists only 
for times t ~ so we can only obtain this estimate on semiclassical time scales. 

This approach is taken in [5, 7, 9, 26-29] and many others, for the wave equations 
and the Schrodinger equations. 

In the proof of the microlocal dispersion estimate Lemma 15. 4[ for the range 
of times t ~ where ^ is localized in a dyadic band ^ ~ 2^, the relation 

(|1.15p implies that (a — [3)/t is bounded by 2^/^, and as a consequence, the kernel 
corresponding to (|1.14[) decays like 2-'/'*t^^/^ on such a time scale. This decay rate 
explains the admissibility condition \1A\ in the main result. The loss in derivative 
comes from taking fjL = 1/4, a = 1/2, and h = 2^^ in Theorem ll.4l 

In Theorem 15.21 we use Theorem 1 1.41 to deduce the semiclassical Strichartz esti- 
mates for (|1.13p as 
(1.16) 

\\U\\ LP ([0,2- 3 /^T])L'i ^ C'(|lC/o||^i/2p -I- ||[/i||^i/2p-,3/2 -I- l|-R|li,i([o,2-3/2T])Hy^''"^^^)' 

where {p, q) satisfies (jl.4p and C > depend on p, q and the Sobolcv norms of V. 



^Here and elsewhere in the paper, the interval [0, 2~^/^T] can be substituted for any interval 
/ of length 2~^/'^T contained in the domain of V, by shifting t to the beginning of the interval /. 

^ To avoid the singularity in the phase at ^ = 0, we assume our initial data arc localized to 
high frequencies. 
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1.3.3. Adding up the dyadic blocks and from linear to nonlinear. We give a brief 
outline of the argument that allows us to move from Theorem 15.21 to Theorem 11.11 
and Theorem 1 1.21 

Let us divide the interval [0, T] into 2^/^ small intervals of the size 2~^/'^T. We 
apply (|1.16p on each short interval of size and we simply sum up 2-'/^ many 

small time-scale estimates. In doing so we introduce an additional loss of l/2p 
derivative. Then by appealing to Littlewood-Paley theory we sum up dyadic fre- 
quencies to assert Corollarv l5.3l We only pause here to remark that the parametrix 
is constructed only for high frequencies; low frequencies can be estimated via energy 
estimates. 

In order to prove the estimate for the nonlinear equation ()1.8p . we employ the 
energy method to establish local existence and uniqueness of the solution of p.Sp 
in the Sobolcv classes. It is detailed in Section [6] Applying 9^ to (jl.Sp . we arrive 
at the linear equation 

5?9> - ^Edldl^u + gUd^d^u + 2udtd^dlu u^dld^u = R{u, dtu) 

for where i? is a collection of lower-order terms. By setting u ~ V{t, a) and 
R{u,dtu) = R{t,a), and applying the above result, we assert that d^u satisfies 
the estimates of Corollary 15.31 As a consequence of uniqueness then u satisfies 
Strichartz estimates as in the Theorem ll.il 

Theorem 11.21 is obtained by repeating the argument above about how to move 
from linear to nonlinear problem via energy method to the Strichartz estimate 
(|1.16p for the dyadic-frequency localization. 

We finally remark that Theorem 1 1 . 21 does not imply Theorem 1 1 . 1 1 since frequency 
localization of the initial data is lost due to the presence of the nonlinearity. 



1.4. Organization. The article consists of three main parts. 

The first part is to formulate the hydrodynamic problem of water waves un- 
der surface tension as a nonlinear dispersive equation. In Section [5] we recall the 
formulation in [3] of the water-wave problem. In Section [3] the system is further 
formulated as a second-order in time nonlinear dispersive equation weakly coupled 
to a transport-type equation. 

The second part concerns the semiclassical Strichartz estimates for the lin- 
earized water-wave equation under surface tension. In Section [J] we construct a 
high-frequency parametrix for each dyadic frequency band and on the frequency- 
dependent time scale. In Section [5] we prove that the parametrix possesses semi- 
classical Strichartz estimates. 

The third part concerns results for the nonlinear problem. In SectionlHl the local- 
in-time existence and uniqueness is established via the energy method. Finally, 
Section [7] presents the proof of the Strichartz estimates for the nonlinear problem. 

Appendix [X] contains a proof of the local smoothing effect for (|1.2p via the 
method of positive commutators, suggested to us by T. Alazard, N. Burq, and C. 
Zuily. Appendices [b1 -ID] collect miscellaneous calculations in the course of the paper 
and linear energy estimates. 
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2. The hydrodynamic problem of surface water-waves 

Recorded here is the approach taken in [3] of the formulation of the water- 
wave problem when surface tension is acted on. The idea is to employ a favorable 
parametrization of the moving surface and choose convenient dependent variables. 

Throughout the paper, partial differentiation is represented either by the symbol 
d or by subscript. The complex plane C is identified with the real two-dimensional 
space E^, whenever it is convenient to do so, via the mapping <I> : K.^ ^ C, y) = 
X ~\- iy. The conjugate of a complex number z is denoted by z. 

2.1. The evolution of the moving surface and the vorticity strength. The 

equation of the moving surface is written as (a;(t, a), a)), where a G M is the 
parametrization of the curve, and <I>(.T(t, a), a)) = z{t,a). Let 

= + Va and 6 = arctan(y„/xa) 

denote, respectively, the square of the arc length and the tangent angle that the 
curve forms with the horizontal direction. The unit tangent and normal vectors of 
the curve are t = (cos 6', sin 6') and n = (— sin0, cos0), respectively. 
The evolution equations of the moving surface are written 

dt{x,y) = [/"t + C/^n. 

In other words, Jj" is the tangential velocity and U-^ is the normal velocity of the 
moving surface. Accordingly, 

respectively. By insisting dtSa ~ 0, and furthermore, = 1 for each t G M_|_ and 
a £ R, we regard the evolution equation of the moving surface as 

(2.1) dte^do^U^ + U^^daO, 

where is determined by solving = U^daO. Such a (renormahzed) ar- 

clength parametrization is assumed initially, and the choice of tangential velocity 
will guarantee that the parametrization is maintained at later time. 

Describing the dynamics on the moving surface, we employ the idea of vortex 
sheets in the two-fluid system, and we suppose that the interface separating the 
vacuum from the fluid moves with different velocities along the tangential direction 
of the interface. 

Let (jy^ represent the velocity potentials of the upper and the lower fluids, re- 
spectively, and let be the densities of the upper and the lower fluids and 
be the corresponding pressures. The Euler equations in the vacuum and the fluid 
region take the form 

(2.2) 9,0± + i|v^±p-KPf = 0, 

2 p± 



The normal velocity U is determined by the equations of motion, while the tangential 
velocity (/H only serves to reparametrize the moving surface. Adding an arbitrary tangential 
velocity does not change the shape of the surface, and thus one may choose the tangential velocity 
to satisfy a certain condition. 
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and the boundary conditions at the interface are written as 

[V0±] • n and [p] = Sd^e, 

where [•] represents the jump of the quantity across the interface. We note that the 
arclength parametrization of the interface offers a particularly succinct expression 
of the mean curvature. 

Let 7 denote the vortex sheet strength. Introducing the Birkhoff-Rott integral! 

(2.3) $(W)(a) = -i.PV r ^^"'^ da', 

27ri z{a) ~ z[a') 

we express the limiting value of velocity at the interface as 

V0^(t, ^^^{z){t, a)) = W{t, a) ± ^j{t, a)i. 

On the other hand, dt{x, y)=W+ (ij" - W • t)t and C/-^ = W • n. 

By combining Bernoulli's equation (|2.2[) with the boundary conditions at the 
interface and by using the above notations, we derive the evolution equation of 7 

(2.4) da = Sdle + 9„((C/II - W • t)7) - 2Wt • t - ^79^7 + 2(C/II - W • t)W, • t. 

The development is detailed in [3, Appendix B]. 

In summary, the water-wave problem consists of (|2.ip and (|2.4p . A useful feature 
of the formulation is that surface tension enters the equation in the linear fashion. 

2.2. The system for the tangent angle and the modified tangent velocity. 

The choice of tangential velocity C/H produces in (j2.4[l nonlinear terms involving 
[/II _ w ■ t. In order to express these terms in a more convenient way, we introduce 
the modified tangential velocity 

(2.5) u=i7-([/ll-W.t), 



and we rewrite the system (|2.ip and (|2.4p in terms of 9 and m, instead of 7. Physi- 
cally interpreted, u measures the difference between the Lagrangian tangential ve- 
locity W-t-|-i7 and tangential velocity J7" which guarantees arclength parametriza- 
tion. Once [x{t,a),y{t,a)) is given, the mapping 7 ^ it is one-to-one. 

The first step is to approximate W in terms of the Hilbert transform. By ex- 
panding <i>(W) in the Taylor fashion, one obtains 

¥(W)(a) =— PV / . 7,1"''' rrda' 



27ri J _^ Za{oL'){a — a') 

+ I 7(a') ( ^T^T^ riV ^ ) 

lmj_^ \z[a)~z{a') Za{a ){a — a ) J 

Note that /C[2;]7 is not singular as the singularities in the expression of JC[z\ cancel. 
Moreover, lC[z\ has the "smoothing" property 

(2.6) WKWWhs ^C{\\e\\Hs+i-r.)\\f\\Hr. for0 1andn = 0,l. 



^ The flow is irrotational. The vorticity, however, has a singular distribution supported on the 
interface. The vortex sheet strength then measures concentration of vorticity along the interface. 

^ In the recovery of the velocity from the vorticity distribution, we employ the Biot-Savart law 
to derive an integral expression, the limit of which at the interface is the Birkhoff-Rott integral. 
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The proof is very similar to that of [2, Lemma 3.5], and hence it is omitted. The 
commutator operator 

TT a -a' 

has a similar smoothing property 

(2.7) mMfW ^C\\h\\H.^,\\f\\H.-s', fors,s'^Oandr> 1/2. 

The proof is found, for instance, in [33, Lemma 2.14]. 
The next step is to represent Wq as 

(2.8) W„ •n= iiJ(7„) + m-n and W„ • t = -iiJ(76l„) + m • t, 
where 



(2.9) 5(m) = =.,t:w|2i-2ipU|i 



zi 



la 



Indeed, by differentiating W = iiJ(7n)+(smooth remainder) and using = —9jt 
one obtains 

Wq = —H{ja)'^ ^ ij^Hl^fOo^i + (smooth remainder). 

The detailed calculation is found in [2, Section 2.2]. 

Using the results above, finally, (|2.ip and (|2.4p are written as 

(2.10a) dtu = ^dle - gO - ud^^u + d-H-r2(t, a)da,9 + (Hda^u + n{t, a)f), 

(2.10b) dtO = -ud^e + Hdau + ri{t,a), 

where 

(2.11) ri(^,a) = -iJ(m- 1) +m • n, 

(2.12) r2{t, a) = Wt • n + uW„ • n + ]p6t + ]pue^. 

The detailed derivation is found in the proof of [3, Proposition 2.1]. 

2.3. Estimates for ri and This subsection concerns the estimates of the re- 
mainder terms in the system (|2.10p . We state the main result. 

Proposition 2.1. The remainders ri and r2 in (j2.1ip and in (|2.12p . respectively, 
satisfy 

(2.13) llnlln. ^Cm\H2A\0\\H^+^){l + \\u\\H^-^) fors^l and 0,1, 

(2.14) llrsllff. s$C(||0||H.+O(l + hllH-^+O' vol. 
Moreover, r2 may be written as r2 = Hdtu + ra, where 

(2.15) llrsllff. 5=: C{\\e\\Hs^^)[l + \\u\\Hs+if for s ^ 1. 

Our result is related to that in [3], but with the important difference that here 
5" > is held fixed whereas in [3] the estimates are uniform as S" — > 0. 
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The remainder term ri involves "smoothing" operators IC[z] and [H,^]. On 
account of (|2.6p and (|2.7p . it fohows that 

(2.16) ||m||ff, s^C(||0||H.+„)ll7lk^- forO 1 andn = 0,l- 
Then, it is immediate that 

(2.17) \\n\\H^ ^Cm\H-^,\\9\\Hs+.)\h\\H2-r. for O 1 and n = 0.1. 

We further estimate ri in terms of u (instead of 7) and 6. Below is the basic 
regularity property of 7. 

Lemma 2.2. Let S > be held fixed. For s ^ 1, if 6 e H^'+i/^ u € and 
7 G then 7 G if" and 

h\\H'^C\\u\\H^+CmHs^U^)- 

Indeed, the definition of u and p.Sp yield that daj = 2daU + H{jda0) — 2m ■ t. 

The assertion then follows from (|2.16p . 

The estimate (|2.13p finally follows by combining (|2.17p with Lemma [221 

A consequence of (|2.16p is that u = ^7 + (lower order terms), which is useful in 

the future consideration. 

Corollary 2.3. For s ^ 1, i/0 G H-'+^^^, u G H' andj e H'-^ then C/H - W-t G 
i/* and 

\\u\\ -w -iWH^, ^cm\H^)\\u\\H.-i + cm\H^)- 

The assertion follows at once from 9q(C/II — W • t) = — • t. 

The estimates for r2 are more involved. Using (|2.8p and (|2.10b ) we write 



r2{t,a) = Wt ■ n + u ( ^-H"(7q) + m • n 

(2.18) 

1 1 

+ -■^{-u9a + HUa + ri {t, a)) + -^udaO. 



Much of our effort to estimate r2 goes to show that the principal part of Wj • n, 
and subsequently, the principal part of 7*2 is Hdtu. 

Lemma 2.4 (Calculation of Wt • fi). For s ^ 1, we have 

(2.19) \\Wt ■ n\\H^, ^ Cm\H^^+^) + C{1 + h||H.+l)^ 

(2.20) Wfn-^H{jt) ^Cm\H^+^) + C{l + \\u\\H^,+.)\ 
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Proof. By writing Wt in terms of 7 and IC[z], we have 
W, -n=Re(!=„«(W,)) 

/I /"^ 7(0;') 1 \ 

+ Re (^-z„(a)PVy_^ ^z.,(a')^^^p^da'j 

- Re ( ^.„(a)PV ( ^) ^44^^da' 



27r \Za{a')J z{a)-z{a') 

■.=^H{jt) +R1 + R2 + R3 + R4. 

In the proof, we examine each Rj, j = 1, 2, 3, 4, separately. 

In order to estimate Ri , we simphfy the expression by introducing 

q3{a, a ) = ; = / Zaa{Ta + (1 - T)a jdr. 



It is standard (see [6], for instance) that 

The Minkowski inequahty and the Fubini theorem then apply to yield that 



00 poo 



7t("') 



2 



9^g3(a, a')\ dada' 



oc J —00 (Q! ) 

whence 

\\R^\\„. ^cm\Hs+^)ht\\L^- 

The smoothing property of IC[z] in (j2.6p implies to yield a similar estimate 

||i?2||H-' SiC{\\9\\Hs + i)H\\L-. 

Further, in Appendix [Bl it is shown that 

(2.21) W^tWHs ^Cm\H^^2) + C\\u\\H^^i for 00. 

By the above estimate for s = 0, then, it follows that 

\\Ri\\h,,\\R2\\h^ ^Cm\H.^,) + C\\u\\H.^i fors^l. 
Next, upon writing W in terms of the Hilbert transform and lC[z], we have 

i?3 = -Re (^Za{a)H f^^atj + iza{a)IC[z] (-^^at 
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whence for s ^ 1 the following inequalities hold: 



\r3\\h^ ^cm\H.^ 



7 



c(||0|U.+O 



-Zat 



L2 



^Cm\H^^^){\h\\HA\Ot\\H^ + h\\m\\Ot\\v^) 

^cm\H.^.){i + \\u\\H.+.f. 

The last inequality uses (|2.10b ). Indeed, 

IMh^ C\\u\\hs+, + Cm\Hs+^) for s ^ 1. 
In order to estimate i?4, similarly, we write 

i?. = Re(i..[i/,..](i-5„(^); 



We claim that 

(2.22) \\zt\\Hs^C\\u\\H^+Cm\Hs+i) fors^l. 

To see this, we write 

zt = $ {(W ■ n)n + (W ■ t)t + (J/ll - W ■ t)t) . 
By (|2.6p and the result of Corollary 12 . 31 then follows 

||W||h. ^^Ch\\Hs+Cm\H^^^)MH^. 

- W-tllff. ^ C{U\H^)\\u\\Hs-.+Cm\H^) 

for s ^ 1, which proves the claim. With (|2.22p immediately follows that 



\Zt\\H' 



da 



Ztda 



L2 



^cmH^^^){i + \\u\\m + \\u\\H^Y. 

Finally, combining estimates for Ri through i?4 yields that 

||i?l + i?2 + i?3 + RA\h' ^ C{\\e\\H^ + ^) + C(l + \\u\\H^ + if. 

This together with (jOTj) asserts ((2T9)) and ((2:20)) . 



□ 



Returning to the estimate of r2, we estimate terms in p.lSp other than Wt • fi 
in the usual way by using (|2.16p . and therefore (|2.14p follows. To establish (|2.15p . 
we write r2 — Hdtu + r^, where 



r3 = dt{U\\ - W • t) + i?i + i?2 + i?3 



Ra 



dtda{U\\ - W • t) = -iff (7,0„) - iff (t^^O + dt{m ■ t). 



Since 



it follows (|2.15p . This completes the proof of Proposition 12. II 

We end this subsection with estimates of the time derivatives of lC[z\f and 
[ff, /i]/, which will be useful in the following section. 



18 



CHRISTIANSON, HUR, AND STAFFILANI 



Corollary 2.5. For s ^ 1 we have 



(2.23) \\dt{IC[z]f)\\Hs ^Cm\Hs^^)il + \\u\\hs+. + ll/llff. + \\dtfh2), 

(2.24) \\dt[H,h]f\\H^ ^\\dth\\H^,\\f\\H^ + \\h\\H^^4dtf\\L^ 



for s ^ 1 . 

The proofs arc in Appendix iBl 

Remark 2.6 (The dispersion relation). We linearize (|2.10p about a flat equilibrium 
u = and 6* = to obtain 



By considering the plane- wave solution u = expik(a — c{k)t), we arrive at the 
dispersion relation 



where c(fc) is the phase velocity corresponding to the wave number k; to = c{k)k is 
the frequency. 

Colloquially, when S* > 0, waves of high frequencies (short waves) propagate 
faster than waves of low frequencies (long waves). Broadening out the wave profile, 
it in consequence induces a certain "smoothing effect". When S" = 0, on the 
other hand, such a smoothing effect is not expected. Alternatively put, the above 
dispersion relation indicates that the linear system of surface water waves exhibits 
a regularizing effect when the effects of surface tension are accounted for. 

The present purpose is to quantitatively analyze such a smoothing effect in terms 
of integrability under the mixed Sobolev norms. 

3. Reformulation: the water-wave problem as a dispersive equation 

The formulation of the water-wave problem under surface tension ultimately 
takes the form of a second-order in time nonlinear dispersive equation, coupled 
with a transport-type equation. 

When both the effects of surface tension and gravity are present, 5 > and 
5 > 0, in (|2.10p the surface-tension term §9^6* is of higher order compared to the 
gravity term gO and it dominates the linear dynamics. For simplicity of exposition, 
thus, the effects of gravity are neglected and further the coefficient of the surface 
tension term is normalized so that 



3.1. Reduction to the dispersive equation. By differentiating (|2.10b ) in the 
i-variable we obtain the second-order in time nonlinear dispersive equation 





q — Q and — 



1. 



d^u - Hdlu = -2udo,dtu - u^dlu - dc,6dlu 

— Zdaudtu — 3u{dau)^ + d^ri + 2udar^ + ur^ + r^. 



Here, ri is defined in ()2.1ip . and 

(3.1) r4 ^d^\r2{t, a)dc,e + {Hd,,u + nit, a)f), 

(3.2) rs ^d-^dt{-H{dtu)d^e + ry,{t,a)d^e + {Hd^u + ri{t,a)f). 
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By ([213]) and ([214]) it follows that 

(3.3) WnWH^ ^Cm\H^+i){l + \\u\\Hsf forOl. 

Further, the leading term of d~^dt{—H{dtu)dad) cancels -dadd^u so that the 
highest-order nonlinear terms in the above equation do not involve 9 explicitly. 
Indeed, successive substitutions of dtu and dt9 by (|2.10p result in that 

r5 = - do.0Hdt{dle - udau + r4) - dtda0H{dtu) 

+ dtrsd^e + r^dtd^e + 2{Hd^u + ri){Hdtd„u + dtn) 

=d^edlu^d~\dl9dlu) 

— idaudtu — 3ud^u + dtTi — udar^ — 2ridau)^ 
+ 9q ^(f3 — Hdtu){Hd^u — daOdaU — udfU — v?daU + dai'i + ur^) 
+ d-\d^Bdtn) + 2d-\Hd^u + ri){Hd^dtu + ^^n). 
Therefore, we arrive at the following equation 

d^u — Hd^u ~ —2udadtu — u^d^u, +R{u, dtu, 9). 
The remainder is given as 

(3.4) R{u,dtu,9) = d'^ri + udaVi + 2ur4 - d^^{d'^u{dtu + udaU - r^)) + rg, 
where 

dafe ~ — {daO)H{—2udadtu — u'^d'^u — daOd'^u ~ Sdaudtu 

— 3ud^u + dtTi — udafA — 2ridau) 

(3.5) 2 2 

+ (r3 — Hdtu){Hd^u — da9daU — udtu — u daU + dari + Mr4) 

+ dc^edtrs + 2{Hda,u + ri){Hdo.dtu + dtvi). 

The remainders ri, r4 and rg involve 6, which incidentally is determined by 
solving (|2.10b ) when u is prescribed. As such, R may be thought of depending u 
and dtu only. In this sense, we write it as R{u, dtu). The remainder term dtu) 
is of lower order compared to udadfU and v?d^u. More precisely, in the following 
subsection we will show that 

(3.6) \\R{u,dtu)\\Hs ^ C{\\u\\H.^^^\^tu\\H^) 

for s ^ 1. 

The water waves under surface tension is finally viewed as the nonlinear disper- 
sive equation 

(3.7) d'^u - Hdlu = -2udadtu - u^dlu + R{u, dtu), 

where R{u, dtu), defined in (|3.4p . is determined with the help of the transport-type 
equation 

(3.8) dtO = -ud^e + Hdau + ri{t,a). 

Nothing is lost in deriving (|3.7p (coupled with p.Sp ) from (|2.10p . To see this, 
we write p.7p as the first-order in time system as 

J dtU + udaU = w, 

I dtv + udaV = Hd^u + R{u, dtu) — dtudaU — ud'^u. 
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Comparing the first equation of the above system with (|2.10|) dictates that the first 
equation of the above is equivalent to (j2.10b .) if we set v = d'^d + r^. It is then 
straightforward to sec that the second equation of the above system is equivalent 
to (|2.10b ) up to constants of integration, which are zero under the assumption that 
the wave profile and its derivatives vanish at infinity. 

Proposition 3.1. The equation (|3.7p . where R{u,dtu) is defined by (|3.4p . is equiv- 
alent to ([2?T0)) . 

Similarly, the initial value problem of (|3.7p prescribed with the initial conditions 
w(0, a) = uo(a) and dtu{0, a) — ui(a) is equivalent to the initial value problem of 
(|2.10p with the initial conditions u{0, a) ~ uo(0) and 9{0, a) = 6*9(0) provided that 
the compatibility condition 

ui = dlOo - uodaUo + d;;^^{~r2{0,a)daOo + {HdaUo + ri(0,Q;))^) 

holds true. 

The transport-type equation (|3.8p is to help to determine certain terms in the 
expression of R in (j3.7p in terms of u and dtu only, and we only need (|3.6p in the 
forthcoming analysis. 

As explained in Section 1, a useful feature of the formulation (|3.7p is that its dis- 
persive character is visible in the linear part. Moreover, the highest-order nonlinear 
terms in p.7p do not involve 9 explicitly. 

Another useful feature of (|3.7p is that it suggests a natural expression for high 
energy of the nonlinear problem. See Section [51 

While (j3.7p is dispersive, its nonlinearity is rather severe, and as such in the 
study of the dispersive property for (j3.7p we must take its nonlinearity into account. 
Indeed, we view (|3.7p as 

d^u - Hdlu + 2udadtu + u^dlu = R{u, dtu). 

That is, we view 2udadtu and u^d^u as "linear" components of the equation, but 
with variable coefficients which happen to depend on the solution itself. Then, we 
make efforts to establish the dispersive property for the linear operator 

- Hdl + 2Vit, a)dtd^ + V^{t, a)dl 

for a general class of functions for V{t, a). 

3.2. Estimates for the remainder. The proof of ()3.6p involves estimates of var- 
ious remainder terms in (|3.4p in terms of u and dtu only (instead of 9). To this 
end, we estimate 9 in terms of u and dtU once and for good. 

Lemma 3.2. For s ^ it follows that 

(3.9) \\9\\hs+2 ^ C{\\u\\h2){1 + \\dtu\\H^ + \\u\\Hs+if. 

The proof is given in Appendix [b1 With the use of (|3.9p . we obtain the estimates 
for the various remainders in terms of u and dtu as 

\\ri\\Hs ^ C{\\u\\H2,\\u\\H^,-.,\\dtu\\Hs-2), 

\\r2\\Hs^C{\\u\\Hs..,\\dtu\\H^), 

WrsllH^ ^Ci\\u\\H.^^,\\^tu\\Hs-l), 
\\r4\\H^ ^ C{\\u\\H2,\\u\\Hs,\\dtu\\Hs-.), 

where s ^ 1. We also estimate for dtri, dtr2, (and in turn, dtr^i). 
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Lemma 3.3. For s ^ 1, 



(3.10) WdtnWns ^ C{\\dtu\\H^,\\u\\Hs^.,\\dtu\\Hs-i), 

(3.11) Wdtv^Wn^ i^C{\\u\\Hs+2,\\dtu\\H^+^), 

(3.12) Wdtr^H^ ^C{\\u\\H^+.,\\dtu\\H^^^), 

(3.13) \\dtr4H^ ^C{\\u\\H^+,,\\dtu\\H^). 



The proofs of (j3.10|) and (|3.12p are given in Appendix [BI Therefore, 

and (|3.6p follows. 

4. Construction of the dyadic frequency parametrix 

This section contains the detailed construction of semiclassical parametrices for 
the linearized water-wave equation. 

4.1. The oscillatory-integral ansatz. Let us invoke the standard notations = 
—idt and Da = —ida and let us denote 

(4.1) P^Dl- iHDl + 2V{t, a)DaDt + V^t, a)Dl, 

where the coefficient function V e iJ'([0, T])iJ'^(]R) is given for some T > fixed 
and for l,k > sufficiently large. We assume that V is real-valued, and we tacitly 
identify V with an HIH^ extension supported in a slightly larger set in t. 
The operator P is obtained by replacing the nonlinear coefficient u in 

d'^t ~ Hdl + 2udadt+u'dl, 

which defines the nonlinear equation p.7p . by a variable coefficient V{t,a)^ and 
thus it is related to the linearized equation of p.7[) . 

Let the frequency cut-off function V'" G C°°(M) satisfy = 1 for |^| ^ A/ + 1 

and i!^{C) = for |^| ^ M for some M > large to be fixed later. Let the dyadic 
frequency cut-off function ij; e C°°(]R) satisfy i^S) eee 1 on f e [2-'^/^ ,2^'^], be 
supported on [2~'^/^, 2^/*], and 

^V(2-^ 1^1)^1 for 1^1 ^ Af + 1, 

where 2^° > M, and let 

^^■(0 = ^(2-^1^1). 

That is, V^XO = 1 on ^ e [2^-1/4, 2^+^'^] and it is supported on [2^-3/4, 2^+3/4]. A 
function / is said to satisfy the dyadic- frequency localization if 

(4.2) i,^{Da)f = f. 

Let j ^ jo be held fixed throughout this section, where 2^" ^ M . Let Mq e L^(]R) 
and u\ e i7^'^/^(R) satisfy the dyadic frequency localization (j4.2p . Our goal is to 
construct a dyadic- frequency parametrix to Pu = with the frequency-localized 
initial data m^'s, n ~ 0, 1, for frequencies comparable to 2^ and on a time scale 
comparable to 2~^/^. More precisely, we shall find a function approximately 
solving 

Pu = in [0,2-J/2T]t X Ma, 
u(0, a) = Uo(q;) and dtu{0,a) — u{{a) 
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(with errors bounded in Sobolcv spaces) for the frequency interval 2^^^ ^ |^| ^ 
2^'^'^. Here, we only consider the time interval t € [0, 2~^/^T], although the results 
apply on any semiclassical time interval of length 2~^/'^T . The over all dependence 
on the coefficient function V is in the fixed time interval t £ [0,T]. For simplicity 
of exposition, we often write |^| ^ 2-' to mean the dyadic frequency band 2-'^^ ^ 
1^1 ^ 2^+2. 

Motivated by the oscillatory integral representation in p.lOp for the zero-coefficient 
case, we make the ansatz 

(4.3) w'{t,a) ^^JJ e-*'^«(e'^"^(*'"'«)/'''+(/3) +e*"''"'*'"'^V'''"(/3))d/3rfe 

for the (leading-order) parametrix. Here, ip^'^ and /■'■^ satisfy the dyadic frequency 
localization (|4.2p . In order to satisfy the initial conditions, we insist that 

^■'"■±(0,a,O 

The phase functions ip^'^ are taken in the class of rough symbols, which is described 
below. 

Symbol classes. For fc ^ and to G M, denoted by 5™ the class of rough symbols 
is defined to be the set 

ST = Ma,0 e {0"'W^^^{R)C^iR) : \d^^ df a\ ^ C,,,™- {0"' for k' ^ k}, 

where (0 = + Symbols in are not necessarily smooth in the a- variable 

(as opposed to classical symbols), but they share in common with classical symbols 
decay properties in the ^-variable. We write Sk for when there is no ambiguity. 

The quantization of a symbol a in 5™ is the usual (left) quantization. It is 
initially defined as an operator on Schwartz functions / as 

Op{a){a,D)f{a)^^ jj a{a,Oe'^"-^'>^ f{(3) d(3 d^, 

and then extended in the distributional sense. We write for the corresponding 
space of quantized operators. 

The main property of the symbol classes 5™ is the boundedness. 

Lemma 4.1. If a ^ S'^ for k > sufficiently large, then Op (a) extends to a 
bounded linear operator 

Op (a) : L^{R) L^{R) 

with the operator norm depending on at most k derivatives ofa{a,S,) measured in 
the L°° norm. 

The proof is exactly the same as in the setting of smooth symbols, keeping track 
of how many derivatives are used. 

Finally, by Wl^'^S'^ we denote the space M^/'°°([0, r])5^, the space of W^'^°°([0, T]) 
functions taking values in rough symbols in 5™. 

The main result of this section is the existence of the dyadic-frequency parametrix 
of the form in (|4.3p with certain properties of the phase functions ip^. 

Proposition 4.2 (Existence of the leading-order parametrix). Let V G -ff'([0, T])H'^{I 
for some T > and /, fc ^ 1 sufficiently large and let j ^ jo ^ 1, where 2^° ^ Al is 
sufficiently large. 
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//(wg,u{) e L^(K) X H ^/^(R) satisfy the dyadic-frequency localization ()4.2I) . 
then there exist the phase functions 

ip^'^it, a, = <± + ^''^{t, a, 0) 

forOs^t^ and 2^-^ < |^| < 2^+^ with i9± e W^2-'T/2t'^°' ^^^'^/i 

satisfy the dyadic- frequency localization ()4.2|1 and f/ie estimate 

(4.4) ||/±|U. «;C7i(h^o|U.+|Iui||^-3/0, 
so t/iat it;^ defined by (|4.3p satisfies 

Pw = E m[0,2-^/^T]t xR^, 
w{0,a) = UQ{a) and dtw{0,a) ~ u-[{a) 

with the pointwise error estimate 

(4.5) Ili^WIU^ < C2{H\\„i + \\u\\\h-i/2), O^t^ 

Here, Ci, C2 > are polynomials in ||l^||//i'([o t])H>'' (K) /"'^ some values of I' and k' 
in the ranges ^ I' ^ I, O^fc'^fc. 

Remark 4.3. As wc will sec later, on the short time seales, the oscillatory integral 
defined by (p^'^ preserves dyadic localization (see Lemma l4.8p . In particular, multi- 
plying the integral by 2^^ , any estimate we prove on w^'^ in has an immediate 
analogue in 

At first glance, the error E looks bad, for it requires in (|4.5p one more derivative 
of the initial data. However, the error is in the inhomogeneity, and as such it only 
appears when we try to measure the difference between the actual solution and 
the parametrix. The energy estimates (j5.4p associated to the linear inhomogeneous 
problem with zero initial data, on the other hand, control 3/2 more derivatives of 
the solution as compared to the inhomogeneity. Hence, E is bad but controllable. 

Remark 4.4 (Remark on the ansatz). We explain why we take (|4.3p as our 
parametrix. 

In order to solve Pu = with the initial data localized to high frequencies, one 
would try a fine parametrix with amplitudes as 

+ e^'p-(t,'^,i)A-(t,a,0f-mdf3d^. 

We require that (f"^ satisfy (^^(0, a,^) = and that A^{0,a,^) and Af{0,a,^) 
are elliptic, and as such the recovery of the initial conditions entails solving for 
a system of elliptic pseudodifferential equations. 

Applying Pu = to our ansatz, we group terms according to their orders in 
^. The worst terms, produced when derivatives fall on the phase functions, make 
a first-order nonlinear equation for (p^, commonly referred to as the eikonal or 
Hamilton- J acohi equation. The other terms form a linear equation, commonly 
referred to as the transport equation, for with coefficients depending on ip^ and 
its derivatives. 

The usual approach to solving the Hamilton- Jacobi equation is through the tech- 
nique of generating functions for the associated Hamiltonian. The equation is, 
however, neither homogeneous nor polyhomogeneous (in ipf' and ^p^), and as such 
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solutions arc found on a time scale comparable to |^|^^^^. See Lemma [4.51 We 
thus construct phase functions (and amplitudes) for each dyadic frequency band 
1^1 ~ 2^ on a time interval comparable to 2"-'/^. 

With (f^ determined, the usual approach to solving the transport equation is to 
expand formal series as 

A±(t,a,0 = E^^'"(*'"'^)r"/' 

and determining recursively. In practice, one takes a finite number of terms 

in the formal series and estimates the resulting error. In our application, we only 
take the very first term, A^'^ = 1, in the amplitude expansion. As a result, we 
have (|4.3p as our leading-order parametrix. 

The full amplitude expansion can be computed assuming more regularity on V, 
but it results in polynomial growth of the lower order terms in the amplitude. 

We now construct the parametrix. Let us fix j ^ jo ^ 1, where 2^° ^ M is 
sufficiently large. Our parametrix as well as the phase function depend on the 
dyadic-frequency band 2^^^ ^ |^| ^ 2^+^. However, we simply write w and f^, 
respectively. To avoid excessive notation, we furthermore consider only the term 
with the superscript -I- and write tp = <^-''~'", / = /-''"'', and 

(4.6) w{t,a) = ^ JJ e-'^^e"^''''"-'^^f{P)df3d^. 

The construction is completely analogous for the term with the superscript — . We 
further assume ^ > 0. After our construction is complete, it will be justified that 
the parametrix preserves the sets 2^^^ ^ C =^ 2-'^^ and —2-'+^ ^ C ^ —2^^^. Thus, 
we tacitly assume ^ > and 2^~^ ^ £, ^ S-'^^ throughout our construction. 
We compute Pw using the ansatz (|4.3p . It is straightforward that 

(4.7) DMt, a) =^JJ e-^<'^e^^{^\ - i^,,)!{P)dl5di, 

(4.8) Dlw(t, a) = ^ // e-^^^e^^Wl " ^^aa)/(/3)d/3de, 

(4.9) DtD^w(t, a) = ^JJ e-'^'^e^'^i^^^t - i^^t)I{Wdt 

We recall that the subscripts mean partial differentiation. In Appendix [Cl wc show 
that iHD^w = —Hd^w can be expressed in a similar fashion as 



HDlw{t,a) 

)f{f3)dpdC 



For simplicity of exposition, in what follows, we take Da and ipa to be positive 
and avoid the absolute value and the sign. Again, this is justified in Lemma l4.8l bv 
showing that negative frequencies and positive frequencies do not interfere. 



STRICHARTZ ESTIMATE FOR THE WATER-WAVE PROBLEM 



25 



Applying Pw = to the results in (|4.7p - (|4.10p . we obtain the following equation 
for ip: 

v'i - Wtt + 2F(i, a){ipa'Pt - iVta - i^Padt - dtda) 

(4.11) 

+ V {t, a){(Pa ~ tfaa) - {^a " ^l^afaa - faaa) = 0. 

In solving the above equation, we only consider terms that are produced when 
first-order derivatives fall on the phase function only. They make a nonlinear equa- 
tion for if, commonly referred to as the eikonal or Hamilton- J acobi equation. For 
other applications, a lower order error may be required. This can be brought about 
by considering a linear transport equation for amplitudes as well; see Remark 14.41 

4.2. Construction of the phase functions. This subsection concerns solving 
the nonlinear Hamilton- J acobi equation and determining the phase functions. 

Lemma 4.5 (The Hamilton- Jacobi equation). Given a coefficient Junction V in 
a bounded subset of H\[0,T])H''{M.) for some T > and for l,k ^ 1 sufficiently 
large, and given j ^ jo with 2^" M > sufficiently large, the following equation 
(4.12) 

ifiUt, a, + 2F(t, a)iptit, a, Ov'aCt, «, + V^t, a)ipl{t, a, ~ vl{t. «, = 
with the initial condition 

has two solutions (p^'^ for 2^^'^ ^ C ^ 2^+^ on the time interval 5j t ^ 2^-'/^T. 
Moreover, 

(4.13) ^''^{t,a,0=ai±e^'{t + ^''^{t,a,0) 

where 

^^'^{t,a,0 = m{\Vit,a)\ + \Vit,a)\')) G 

satisfies 

\dfd^^d^\ < C™'2--'/2 (0"™' fork' ^k~l. 
The derivatives of f^'^ have the following properties 
(4.14) 

(4.15) 

ip^^{t,a,i)^±e^\l + §f(t,a,i))^e^\l + 0(C"')) e W^2-T/.T'5'-i 

for 2^~^ ^ ^ =^ 2^+^. That is, (p^^^ is of order ^ on the dyadic frequency band 
2^^^ ^ C ^ 2-'+^ and on the small frequency- dependent time scale 5j i ^ 2^^^^T 
and (fl'^ is of order ^^/^ . 

The class of rough symbols 5™ is defined in the previous subsection. 
We recall that only (p = (y9^'+ is considered in writing (|4.1ip . That enforces ^ be 
positive, which is tacitly assumed in this subsection and the following one. 

Remark 4.6. A few comments are needed about Lemma l475l 

First, the result applies equally well on any time interval / of length 2~^/'^T 
which is in the domain of V . The initial conditions for ip^ are then prescribed at 
the beginning of / rather than at 0. 
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Second, since ()4.12|) is quadratic in Lpt, there are two phases ip^. This may 
be thought of as an analog of incoming/outgoing solutions of the wave equation, 
although we do not exercise that level of sophistication here. We write the two 
branches of (|4.12p as 

(4.16) ipf^~V{t,a)vt±{vt.f'\ 

and we will work on the "factorized" Hamilton-Jacobi equation in sequel. 

A standard method of existence for a Hamilton-Jacobi equation of the kind of 
(|4.16p would be to construct Lp^ as a generating function of a symplectic transforma- 
tion which arises as a solution of the corresponding system of ordinary differential 
equations 

fd = -y(i,a) + |r;i/2, 

with the initial condition 

a(0) = /3, ry(O) = C 

for each /?, C £ M. Here and elsewhere, the dot above a variable denotes the differ- 
entiation with respect to the i-variable. This system has a solution only for a time 
scale comparable to which is what we are after, but we desire a finer control. 

The idea lies in that (|4.16p is sought of as perturbation of Lpf = (y'^)'^^^ under 
a lower-order term. With the initial condition (p(0, a,^) = the solutions of 
this equation arc found explicitly to be iy9o(t, a,^) = ± t^'^^^. Then, it seems 
reasonable to find solutions to (j4.16[) as a perturbation of Lp^. The added term in 
()4.16p . while being of lower order, destroys the homogeneity of the equation and it 
causes serious difficulties in the application of the Hamilton-Jacobi theory. 

Proof. For simplicity of exposition, we will suppress the dependence of the phase 
on dyadic frequencies f ^ 2^ and we will prove for the + sign only. Let ^p = (^^'+ 
denote the solution; the proof for the — sign is identical. 
Let us consider the initial value problem 

(4.17) ^t = -V{t,a)^^ + ip''J\ v'(0,a,O=aC, 

where 1^9 is a function of t, a with a parameter ^. As is remarked above, we construct 
the solution as a perturbation of 

(^o(i,a,C) + 

3/2 . 

which is a solution to a homogeneous equation Lpt = (pj with the same initial 
condition. Specifically, we make the ansatz 

ipit,a,0=aC + e^^{t + dit,a,0). 

where 

i?(i,a,C) = i5(i,2-^/^a,0 and ^ e Sk- 
Substituting our ansatz for ip into (j4.17p yields the following initial value problem 

.4 = -r'/'Vit, 2^/-'a){l + + (1 + 2-^/2^1/2^^)3/2 _ 1^ 

\^(0,a,0 =0 

where "d is evaluated at (t,a,^). The corresponding Hamiltonian is 

q(t, a, rj) - -r''^V{t, 2^/^a){l + 2-^'^e'^v) + (1 + 2--'"/2^i/2,^)3/2 _ ^ 
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and the corresponding Hamihonian system is 
(4.19a) <^ 

with the initial conditions 

(4.19b) a(0) = (3 and 77(0) = C, 

where G M and C G for some e > 0. We recall that the dot above a 

variable denotes the differentiation with respect to the t-variable. Since under the 
assumption of ^ ~ 2^ the right sides of (|4.19b .) satisfy the Lipschitz condition with 
respect to a and 77 with the Lipschitz constants comparable to 2-'/2, it is standard 
from the theory of ordinary differential equations that a unique solution of (|4.19p 
exists on the time interval ^ i ^ 2~^/'^T for the range of the initial conditions 
given above, and the solution is at least as smooth as the right hand side. 

We write a = a*{f3, C)j = V*{0i C)j and «*(/?, Q = {a, 77) for the symplectomor- 
phism given by the solution of (|4.19p . Let uj be the 1-form 

ui = —rdt + rjda + (3dC,, 

and let A be the surface 

A = {(t, q(t, 0), C), : /3, C) G [0, 2-^/^T]t x x [-e, e]^}. 

Since A is a graph it is an embedded three-dimensional submanifold of r*M''. Since 
dw is a symplectic structure on T*M?' and since the fact that k* is symplcctic 
implies dw|A = 0, additionally, A is a Lagrangian submanifold. Such an embedded 
Lagrangian submanifold A can be written as the graph of a closed 1-form, say 
a{t,f3,()- Since M.^ is simply connected, the Poincarc lemma implies that there 
exists i?(t, /3, C) such that 

d^ = a. 

We will prove in Appendix [Cl that the mapping 
(4.20) P^a*iPX) 

is invertible for each ( G [— e, e] and ^ i ^ 2~^^^T, and we write /? = P{a,Q. 
Comparing uj with dd written in the (i, a, C) coordinates, we find that 

dd = -rdt + rjda + /3dC- 

This implies 

dd dd 

— =/3(a,C), = ?7(/3(a,C),C), and —^q{t,a,r]), 

oQ oa at 

which, in turn, implies for C G [— e,e], 

Mt: «, C) = -C'^^V{t, V'^a){l + 2-^-/2^1/2^^ 

+ (l + 2-^/Y/'^o)'/'(i,a,77)-l. 

Therefore, (|4.17p follows once we substitute for 1? and ip. 

The estimates (|4.14p and (|4.15p follow easily by differentiating the relations 
(|iT7)l and dHH), using the timescale ^ t < 2-^/'^T and substituting i?(i, a,C) = 
^{t, 2-J/2a, 0- This completes the proof. □ 
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4.3. Finishing up the construction: recovery of initial conditions. It re- 
mains to determine /^'* to satisfy the initial conditions 

w-'{0,a) = Uq(q;), Dtw-'{0,a) = — m{(a). 

Since (p^'^{0, a, ^) = it follows that 

A\-+{a, D^)P-+ + {a, D^)p--. 
Thus, we are led to the elliptic system for f^'^ as 
P-+ + P- = ul 

A{'+{a,D^)P'+ +A{'-ia,D^)P'- = -iu{. 
The results of Lemma 14.51 state that 

behave like the classical symbol ±^'^/^. Consequently, A{'^ are elliptic pseduod- 

3/2 

ifferential operators, which are approximately ±Da , and they have approximate 
inverses, denoted by {A{'^)^^ . We only pause here to note that the existence of 
approximate inverses here means the existence of honest inverses with the same es- 
timates as the approximate inverse maps. Indeed, the error involved in our setting 
is 0(^~^/^), which, in the dyadic-frequency band ^ ^ 2^ is 0(2^-'/^). Hence the 
approximate inverse mapping has an error which is bounded by 0(2"''/^) on the 
appropriate Hilbert space and the honest inverse can be obtained by the natural 
Neumann series. 

Therefore, the operators 1 — {A-['~p^A\''^ and 1 — {A-'{'^)^^ A\"~ are bounded on 
and invertible with bounded inverses. We set 

p+ =(1 - (Ar)-M{'+)-iK + *(^r)"'"i), 

P^- =(1 - (Ai'+)-Mr)-^K + ^iA\:+r'ui). 

Then, l|4.4p follows immediately. 

This completes the construction of the paramctrix and the proof of Proposition 
14:21 

4.4. The Fourier integral operators. We establish the basic boundedncss prop- 
erty and propagation of singularities for our leading-order paramctrix. Let us first 
establish a basic L^-boundedness result. 

Lemma 4.7. Let F{t) for $J t ^ 2^^/^T be defined, initially on the Schwartz 
class, by the formula 

F{t)f{a) = JJe-^^^e''^'^'"'^^fiP)df3d^ 

where </?-' is either of ip^'^ and is constructed in Lemma \4-5\ Then F{t) extends to 
a bounded linear operator L\ and 

\\F{t)\\Li-.Li 1 + 0(0- 
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The proof of the lemma shows that F{t) is almost unitarity point-wise in t. If 
(ys^ have been constructed on an interval of length 2^^^^T in the domain of V{t, a), 
then the result holds equally true for each point on the interval. 

Proof. We will prove the assertion for ip' = </5-''^, the proof for (^■''~ being analogous. 
Further, We will prove the estimate for F* rather than F, since some nontrivial 
cancellation makes the proof much easier. Indeed, we will show that F*{t) ~ 
F~^A{a,2~^/'^Da) for a pseudodifferential operator A with rough symbol, A = 
1 + 0{t). Since 

= {F{t)F*{t)fJ) < \\F{t)F*it)fh.\\fU.^, 
it suffices to prove the assertion for 

(4.21) F{t)F*{t)f{a) ^ JJ e*('^'(*'"'«)-'^'(*'"'^«»/(a')rfa'de. 

Recalling the results of ip from Lemma 14.51 we have 

^^{t, a, - V'{t, = {a- + {t, a, - ^\t, a', £)) 

= (a-a')C(l + ^''(i,a,a',C)) 
where §^ = &{t) e w!^'^/^^Sl_^ satisfies 

for ki + k2 — 1. We perform the change of variables rj = ^(1 + a, a', C)) in 
to obtain 

F{t)F*{t)f{a)= 11 e*("-"')«A(t,a,a',0/(a')d«'rf^, 
where a symbol 

satisfies 

d^'d^^A = 2-^('=i+'=^'/20(^-) I <^ ki + k2 k - 1. 

Then the Calderon-Vaillancourt theorem implies the assertion. □ 

The following lemma regarding the Fourier integral operator related to (|4.3[) 
shows how to pass derivatives through the oscillatory integral, and it also justifies 
localizing to positive or negative ^. 

Lemma 4.8. Let to S M, / ^ 2, fc ^ 3 and let j ^ jo ^ 1 with 2^° ^ M sufficiently 
large. Let ip^'^{t,a,^) € '^'^ constructed in Lemma \4^ 

Suppose that B{a,^) £ ^-^^ where k' ^ k + 4. Then, for any f e i?™(K) 
satisfying the dyadic frequency localization (j4.2p the following 

Bia,D^) yye-*'5«e'^"*(*^"'«)/(/3)d/3de 

-'Ke"^''^'-*^°''^h^"'B{t, a, e)/(/3) d(3d^ + {Ef){t, a) 



holds, where B G W^j-jV^r'^fc-i supported in Cg2-' < ^ ^ Cj2-' for some < Cq < c'^ 
independent of j, and Ef satisfies the dyadic-frequency localization (|4.2p and Ef 
enjoys the estimate 

II-^/IIl2([o,2-j/2T])L2 < Cl|/l|^™-i. 
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Furthermore, ifip^{Da) G ^'°/ is equal to 1 in the dyadic region 2^^^ ^ |^| ^ 
2J+^ then 

i^HD^) J J e-'^^e'^"^^'''''^^fiP)dpd^^ JJ e-"3«e'^"*(*'"'«)/(/3)d/3de 
modulo a lower order error. 

As usual, the constant C > in the error estimate depends on up to k derivatives 
in a of B{a,£,) and k derivatives in a of V(t,a). The requirement of 4 derivative 
comes from keeping track of the number of derivatives used to control the error 
terms in the Egorov theorem. Of course, the error estimate can be improved upon 
more careful use of the Egorov theorem. 

Proof. As usual, wc prove only for ip^ = the proof for ip^'~ is similar. 
Let us consider the oscillatory integral operator 



The results of Lemma 14.51 say that 

^^{t,a,0^< + e^'{t + i)'{t,a,0) 
with d^{t,a,^) = 0(t) G M^^L^/^T^fc and that 

Since 



it follows that 



(4.22) r s; C 



dip- 



da 



2 



dip- 



2^ 



for 2^-2 <; ^ ^ 2^+2, 



where C > is independent of ^. In light of this and Lemma 14.71 then F is an 
elliptic Fourier integral operator. 

Let us also choose X2(0 such that X2{£,) = 1 on [1/2,2] and it is supported in 
[1/4,4]. Let xiiO = X2(2~-'0- We define a modified phase function 

(4.23) ^{t,a,0 = a^ + te^' + xm^'it,c.:m'^'- 

Note that the modified phase function (p^ is defined for all ^ ^ M and (f^ {t, a, ^) = 
(p^{t,a,^) on the dyadic- frequncy band 2^^^ =^ C =^ 2-'+^. 

We note that the phase ip^ is the generating function of the symplcctomorphism 
in the proof of Lemma [4.51 in the dyadic band 2^~^ ^ C ^ 2''+^, which is a lower 
order perturbation of the symplcctomorphism 

(4.24) a^->a+^t^i/2, 

The phase ip^ generates the same symplcctomorphism in the dyadic region, and 
extends it to be (|4.24[) in the rest of phase space. Let k* be this extended symplcc- 
tomorphism. 

In light of the Egorov theorem [13], then it follows that F transforms symbols 
according to the symplectic transformation . 
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It remains to show that maps dyadic frequencies to dyadic frequencies and 
preserves the order of the symbol. Indeed, the ^ component of k* is ^(1 + 0(t)), 
whence 

{co2^ sC e ^ ci2^} C {{^%{a,0 ■■ c'^V ^ £, < ^2^} C {c'^V ^ ^ < c'lV} 

for some positive constants cq < cijCg < c'j^, and Cq < c", where (k*)2 denotes the 
second component of k*. 

Therefore, for any pseudodifferential operator i? G vj/™^ where k' ^ fc + 4, it 
follows that 

B{a,D^) J J e-'^^e"^'^*'"'^^f{Odpd^ = JJ e-*''«e''^' (B/)(/3)dM 

for some pseudodifferential operator B G with principal symbol 

a(B) = e(i,a,0(^*)V(i?), 

where e £ S^-i is elliptic on the support of {hi*')*a{B). This completes the proof. 

□ 



Remark 4.9. The assertion of Lemma 14.81 holds true when replacing ^ by — ^. 
This justifies considering only positive ^ in the construction of the parametrix in 
the previous section. In what follows, we will drop the assumption that we work 
on positive ^. In other words, the phase functions take the form 

and similarly for the amplitudes and others. 



5. StRICHARTZ estimates for the LINEARIZED EQUATION 

A dispersion estimate on a semiclassical time scale implies semiclassical Strichartz 
estimates for linearized water-wave problems under surface tension. Gluing these 
together Strichartz estimates on a fixed time scale are established with loss in de- 
rivative. 



5.1. Preparation for the proof. Let us consider the initial value problem of the 
linear homogeneous equation 



(5.1) 



dju - Hdlu + 2V{t, a)dcdtu + V'^{t, a)dlu = 0, 
M(0,a) = Mo(a) and dtu{0,a) = wi(a), 



where the coefficient function V G iJ'([0, T])i?^(M) is given for some T > fixed 
and for /, fc ^ 1 sufficiently large. Let us also consider the related initial value 
problem of the inhomogencous equation with the zero data 



(5.2) 



d^v - Hdlv + 2V{t, a)dcdtv + V^{t, a)dlv = R{t, a), 
v{Q,a)~Q and 9tu(0, q;) = 0, 



where R G ^^([0, T])iJ^(M) for some s ^ 0. 

As a preliminary result, we establish the existence and uniqueness of the actual 
solutions of (|5.ip and (|5.2p via the standard energy method. 
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Theorem 5.1 (Existence and uniqueness for ([51]) and ^^). Let V G i/'([0, r])iJ^(M) 
for some T > and for some l,k > 0. 

For each pair of uq G H''+^/'^{R) and ui G H^iR), where 0^s + 3/2^fc, there 
exists a unique solution u to (|5.ip on the interval ^ t ^ T satisfying 

(5-3) ll"llicc([o,T])ii-;+^/^(R) + ll'^*"lli-°°{[0,T])HJ(R) ^ C'i(II"o|Ih= + 3/2(r) + ||ui||^s(k)). 

Furthermore, for each R G L'^([0, T])_ff^(R), there exists a unique solution v to 
(|5.2p satisfying 

(5-4) lkt||L°=([0,T])ffJ(K) + lkllL=c([o,T])ff;+^/=(B) ^ || l2([o,t])//j (R) ■ 

ffere, Ci,C2 > are polynomial in \\V\\L^^[o,T])wi-°°{R} arid \\Vt\\L^(^[o^T])wi'°°{R)- 

The existence and uniqueness is standard by combining the energy estimate (|5.3p 
and (j5.4p with rcgularization and a Galcrkin approximation. The detailed proofs 
of ([5?3|) and arc in Appendix [d1 

The main result is in this section is scmiclassical Strichartz estimates for the 
linearized water-wave problem under surface tension. 

Theorem 5.2. Let V G i/'([0, r])i?^(R) for some T > and for /, fc > 1 suf- 
ficiently large and let j ^ jo, where jo 3> 1 is sufficiently large. Suppose that 
U G L°°([0, 2^J/2r])iji/2p(K) solves 

(5.5) 

d^U - HdlU + 2V{t,a)dadtU + V^{t,a)dlU ^ R{t,a) for t ^ [Q,2-'j/^T], 
U{Q,a)^Uo{a) and dtU{0, a) ^ Ui{a), 

where {Uo,Ui) G H^'^p{R) x H^/'^p-^/'^{R) and R G L'^{[0,T])H^/'^p-^'^{R). Sup- 
pose, further, that U , R, Uq, and Ui satisfy a dyadic frequency localization (|4.2p at 
frequency 2^ . Then, U enjoys the estimate 
(5.6) 

l|f^llLP([0,2-J/2T])L« ^ C(l|t^o||^i/2p + ||C/l||^i/2p-3/2 + l|-R|lii([o_2-J/2T])ify^''^^''^)' 

where (p, q) satisfies 

2 1 _ 1 

p + g" 2' 

and C depend on p, q and the Soholev norms of V of orders bounded by I and k. 

Moreover, gluing together these estimates on 2^/'^ many time intervals of length 
2-j/2y TffQ obtain Strichartz estimates on a fixed time scale with loss in derivative. 

Corollary 5.3. Let V G H\[Q,T])H^R) for some T > and for some /, fc > 1 
sufficiently large. 

If{Uo,Ui) G H^/P{R)xH-^/^+^/P{R), then the solution U G L°°([0, T])iJi/P(M) 
of the initial value problem (j5.ip satisfies the inequality 

(5-7) \\U\\lp{[q.t])l^{r) =^ C'i(l|t^o||ifi/p(R) + ll'^i||_H--3/2+i/p(R))- 

IfR G L2([o,T])iJ-3/2+i/p(M)^ then the solution v G L°°([0, T])i?i/p(M) of the 
initial value problem (|5.2p satisfies the inequality 

(5-8) lkllLP([0,T])L<!(R) ^ C2||-R1Ili([o,T])//-3/2+i/p(R)- 

Here, {p,q) satisfies the admissibility condition in Theorem \ 5.2\ with q < oo, and 
Ci, C'2 depend on p, q and the Sobolev norms of V of orders bounded by I and k. 
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5.2. Proof of Theorem 15.21 and Corollary 15.31 We explain our strategy to 
proving Thcorcm l5.2l We first demonstrate that the homogeneous dyadic-frequency 
paramctrix in Proposition 14. 21 satisfies an L]^ — > dispersion estimate on a semi- 
classical time interval of length 2^-'/^. Then we use Theorem 11.41 to deduce the 
semiclassical Strichartz estimates for the parametrix. We use the energy estimates 
to show that the parametrix is sufficient to estimate the actual solution in the 
Strichartz norms on semiclassical time scales. To prove Corollarv 15.31 we sum up 
the estimates on semiclassical time scales to obtain an estimate on a fixed time 
scale with loss in derivative from the summation. 

We begin by proving the short-time dispersion estimate. Without loss in gen- 
erality, we work on the time interval t G [0, 2~-'/^r]. The analogous statments 
are true for any interval of length in [0,T]. In order to estimate various 

quantities involving the parametrix (14. 3p we have found it expedient to "glue" to- 
gether on the time scale ^ t ^ 2^^^^T the phases (p-''^, which are supported on 
2"'"^ ^ ICI ^ 2-'+^, to obtain global phases with similar properties. As usual, we 
avoid excessive notation by considering the (p^'"*" case only and writing tp' = 
and similarly for other quantities with the ± parity. The analysis in the (p^'~ case 
is completely analogous. 

For e > small, let us choose Xi(^) such that xi(t) ~ 1 on [0, T] and it is 
supported in [0 — e,T + e]. Recall that X2(C) is such that X2{0 = 1 on [1/2,2] 
and it is supported in [1/4,4]. We need X2 satisfying x^iO = 1 on suppx2(C) 
and SUPPX2 C [l/4-e,4 + e]. Let xiW = Xi{2^^H), xHO = X2(2--'0, and 
xiiO — X2(2~-'0. We define a modified phase function 

(5.9) ^{t,a,o = a^+m^^^+x{{t)xm)^'it,a^m^\'^', 

where I?-' ~ has been constructed in Proposition 14.21 The modified phase 
function ip^ is defined for all t and all |^| ^ A/ and we claim it has the properties 

(5.10) ^'{t,a,0 = torn + ler^/^) G W^L%^S°. 
To see this, note that 

(5.11) m.c^,o = o{\t\ + \(r'/'). 

Indeed, a simple substitution of ip in the eikonal equation yields that satisfies 

(5.12) Mt, a, = -\^r'^^V{t, a){l + + (1 + d^f/^ - 1, 

where i}a{t,a,^) = 0{t). Since i?(0,a,O = by construction, i?(i,a,0 = ^0(1^1 + 
^1^^/^) as claimed. 

We redefine w^'^ by replacing the phase function ip^ by ip^ . The new oscillatory 
integral agrees with the old one for ^ t ^ 2~-'/^T and for the data Uq and u\ 
satisfying the dyadic- frequency localization (|4.2p . but it has the virtue of being 
globally defined so that we can use the theory of Fourier integral operators. 

Let us consider the Fourier integral operator 

(5.13) {Fm,a) ^^11 e-'^«e'^^^(*'"'«)xi(t)x^2(0/(/3)dM- 
It is standard from the results of Lemma H771 and Lemma that 



IIW)lli/i<(i + 0W)ll/ll//s 
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for s ^ and each ^ t ^ T. As in Section [L3l for the Schrodinger equation, we 
write the oscillatory integral as 

Ff{t,a) = J K{t,a,[3)f{(i)df3, 

where 

(5.14) K{t,a,(3)= J e-'^^^e^^'^*'"'ihimmdC 

Lemma 5.4 (Microfocal dispersion estimates). The kernel in ()5.14p satisfies 

\\k\\l^ C2^'/4t-i/2 for \t\ s: 2-^/2r, 

so that 

\\Ffit)h^ ^ C2^/H-^/^fhi^ for \t\ ^ 2-^/2r. 

Proof. The proof uses the method of stationary phase together with the assumption 
that t is in an interval of length to estimate the error terms. We assume as 

usual that ^ is large and positive; the large negative case is treated similarly. 
We next compute 

+x{{t)xume^'d^^={t,a,o 
^a~(3+^t\e/'+t\e^'0(\t\ + \^\-^/') 

on \t\ 2-^/2y foj. 2J-1 <^ ^ 2^+\ The second equality uses that e 5^, 
the t localization due to Xii^), and (|5.10p - (|5.1ip . A critical point of the phase 
(f^ {^;t^ a) — (3^, therefore, is given by 

We observe that the localization of X2 implies that the integrand is zero unless 
^ ~ 2-' , and hence ~ 2^ ■ 
Next, we compute for \t\ 

di{^^{e,t,a)^(3o = lm-'/' + ^xmxmm-'^'^'it:c.:0 

+x{{t)xmm'^'dl^'it,c.,o 

+ 3xi(i)2-^X^(2-^|e|)|^r/2z9-'"(i,a,0 
+ xl(i)2--'x2(2-^|CI)l^l'/'a^i?^(i,a,C) 

= ^t\c\-'/'{i+0{\t\ + \cr/^)), 
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where we have again used the t locahzation of xi(0 ^^^^ (|5.10|) - (|5.11|) . Evaluating 
this at the critical point, we obtain 

8 p — a 

Thus, the critical point is nondegenerate for i > 0. 

Applying the method of stationary phas^ completes the proof of the lemma, 
once we observe that the restriction \t\ ^ 2~^^'^T combined with the restriction 

~ 2^ implies that (/3 - ay/^/t^^^ is bounded by 2^/^ for this range of t. □ 

In order to apply Theorem 11.41 we need to control F{t')F* (t). 
Lemma 5.5. The Fourier integral operator F given in (|5.13p satisfies the estimate 
\\F{t')F*{t)f{a)h^ ^C2^/^\t' -t\-'/^f\\L.^ 

for \t'\, \t\ < 2-i/^T. 
Proof. We write 

F{t')F*{t)f{a) = |y'e*('^'(*''"^«)-^'(*'"''«»Xi(t)xi(t')(xKO)'/(«')rf«'rfe 
and its corresponding kernel 

We cannot apply the same argument as in Lemma 15.41 to estimate Ki , since the 
difference of the two phases depends on a and a'. However, the only part from 
which we cannot directly factor the t' — t is of lower order. Precisely, 

^\t',a,0-^\t,a',0^ia-(^')^Ht' -mi'^' 

= (a - + 0(2"^/')) + (f - + 0(2-^/')). 

We then apply the same arguments as in Lemma 15.41 to complete the proof. □ 

The parametrix constructed in Proposition |L2] is not an exact solution for (|5.ip . 
and it only exists for a short time scale depending on the dyadic frequency band. 
We use energy estimates for the linear equation (|5.ip to show that the parametrix 
is sufficient to estimate the actual solution. 

Proof of Theorem \5.2l We break U into a piece solving the homogeneous equation 
with the frequency-localized initial data and a piece solving the inhomogeneous 
equation with the null data. Once we prove the result for the homogeneous equation, 
then the result for the inhomogeneous equation follows immediately from Duhamel's 
principle. Hence, we may assume that U solves the initial value problem for the 
homogeneous equation 

PU^o for t e [o,2-J/2r], 

f/(0,a) = [/o(a), dtU{0,a) = Uiia) 

^Stationary phase gives a bound of 2^^'^t~^^^ plus a reminder term. Since we localize both in 
time and space this reminder is even smaller than t~^^^. 
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recalling the notation (j4.ip . We approximate U on the interval [0, 2^-'/^r] by the 
oscillatory integral parametrix given in Proposition [421 The parametrix W solves 

J PW = E for t e [0, , 
\w{0,a) = Uo{a), dtW{0, a) = Ui{a), 

where E satisfies the error estimates (|4.5|) . We write 

\\U\\lp([0,2-^/^T])LI ^ II^IIlp([0,2-J/2T])L« + " ^IIlp([0,2-j/2T])LJ ■ 

Upon applying Theorem 11.41 it is readily seen that the first term on the right side 



satisfies (|5.6p . 

The function Z = U — W solves 

jpZ = -E for t e [0, 2-^/'^T], 
\z{0,a) = 0, dtZ{0,a)=0, 

and hence satisfies the energy estimate (|5.4p 

(5.15) ||9tZ||ioo([o^2-3/2T])_f/5 + \\^\\L'=°{[0,2-i/^T])Hi+^^^ ^ (^\\E!\\l^{Io,2-i/2T])H^ 

for -3/2 s < fc-3/2. 

By the Sobolev embedding, Holder inequality in time and the energy estimate 
with s = 1/2 — l/(j we then have 

II^IIlp([0,2-3/2T])LS LP([0,2-J/2T])//y^-i/« 

^C(T2-^/2)i/p||£;||^,^j^^^_^^,^j^^_,_,,, 

Then using the estimates on E given in (|4.5p . the estimate (|5.6p for Z follows. □ 

Proof of Corollary ] 5. S[ As above, once we prove (|5.7p . then an application of Duhamel's 
principle and the Minkowski inequality yields (|5.8p . In what follows, thus, we con- 
sider only the homogeneous equation. 

The first step is to chop the actual solution of (|5.ip into pieces each of which is 
localized in a dyadic frequency band. Let us choose a partition of unity in ^ 

i = (i-v.°)(o + E^'(^) 

as in Section m It is standard from the Littlewood-Paley theory that if / G L''(M), 
q < oo, then 

Ii/IL, ^ cii(i - i^'){D^)fhi + ( E ' 



and hence it suffices to prove (|1.3p on dyadic frequency bands. 

Let U be the actual solution of (|5.ip . and let = ip^{Da)U . It is readily seen 
that solves 



(5.16) 



d'^U^ - Hd^W + 2V{t, a)dadtW + V^{t, a)dlW R^iU), 
W{0, a) = U^{a) and dtW{0, a) = U({a), 
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?A^pa)[/o, Ui = ^'J{D^)Ui and 

W{U) = [2V{t,a)dc.dt + V^{t,a)dl , i''{Do,)]U 

= 4''{Da)2-^A2v{t, a, DMDa + Ay^ {t, a, D^)Dl)U. 

Here, [•, •] denotes the commutator, A2V and are zeroth-order pseudodifTeren- 
tial operators, and ijj^ is a smooth function with support contained in a neighbour- 
hood of the support of V'-' ■ It is immediate to see that is comparable to 2^ on 
the support of tp^ . Consequently, we have the estimate 

(5.18) < C{UW\\hi^ + U^dtUWLi), 

where C > is a constant independent of j depending only on a finite number of 
derivatives of V . 

Let [/-' = Ul + Uf be the solutions to the homogeneous and inhomogeneous 
problems, respectively, corresponding to (|5.16p . That is, Uj^ solves 

PUl = for te[Q,T], 
Ui{0,a) = U^{a), dtUliO,a) - UHa), 

while [// solves 

PUl=m{U) for tG[0,r], 
f//(0,a) = 0, 9tC/f(0,a) =0. 

We will prove the Strichartz estimates for [/^, which, via Duhamel's principle, will 
imply the Strichartz estimates for Uf . Then summing in j will imply the Strichartz 
estimates for U solving (|5.ip . 

We divide the interval [0, T] into 2^/^ small intervals of the size 2~^^^T. Let 
j = 2-il'^(m, - 1)T, where 1 < m 2^/^ and let I™-^ = [T,„j , r,„+i,j] so that 

[0,T]= y 

l^m^2i/2 

Then, we apply the results of Theorem 15.21 on each short time interval 7™'-', 1 ^ 
< 2-'/2, to obtain 

ll^/illLP([0,T])L5!(E) ~ 5Z II^/iIIlp(/'".J)-L9(R) 

771—1 

^ C ^ (l|C^/l(T'77l,i)llifl/2p(R) + ||9tt//j(T„ij)||^l/2p-3/2(R)^ 

m— 1 

_ffl/2p-3/2(R) I . 

The last inequality uses the energy estimate. 

On the other hand, 2-'/^^ ~ d]J'^^ on the frequency support of [/q and J/f, and a 
localized version for j > jo of (|5.7p follows. To finish the proof we first notice that 
for small frequencies, j ^ jo, the estimate follows from Sobolev embedding and the 
energy estimate (j5.3p . By Duhamel's principle and Minkowski's integral inequality 
then also (|5.8p follows. This completes the proof. □ 



where ~ 



(5.17) 



38 



CHRISTIANSON, HUR, AND STAFFILANI 



6. Local well-posedness via energy estimates 

This section concerns the local well-posedness of the initial value problem asso- 
ciated to (|3.7p . The result is of independent interest. 

We recall that R{u,dtu) (|3.7[) . defined in (|3.4p . involves 9, which is determined 
by ()3.8p . To be precise, hence, the local well-posedness is to be established for ()3.7p 
coupled with (|3.8p . Upon the observation that ()3.7p and ()3.8p are of different type, 
we proceed by a "bootstrapping" argument. First, given 6 in an appropriate func- 
tion space, the local well-posedness for p.7p is established via the energy method, 
with R{u, dtu) evaluated with the given 9. With u so obtained , next (|3.8p is solved 
via the standard method of characteristics with the variable-coefficient u and ri 
evaluated by the solution u in the first step. We focus on the development of an 
energy estimate for ()3.7p and its local well-posedness. 

We begin by writing (|3.7p as the first-order in time system 



(6.1) 



J dtU = V — udaU, 
]dtv = Hd^u — udaV + R{u,v). 

In other words, v = dtU + udaU is the material derivative of u. Here, 

R{u, v) — R{u, V — udau) -f vdaU + uidau)^ 

satisfies an estimate similar to p.6p : 

\\R[u,v)\\hs ^C[\\u\\Hs^^^\v\\Hs) forO 1. 

Let us define the fc-th energy associated to (|6.ip as 

(6.2) £'^{t) = ^ / {{d^+'u)Hd^{d^^+^u) + (a^«)2) da 
and the energy function for (|6.ip of order s as 

s 

(6.3) r(t) = h||i.(t) + ||z;||i.(t)+^f'^(t). 

fc=i 

Note that the operator Hda is a positive operator with the symbol of its Fourier 
transform |^| and that 

/oo 
{f + fHd^f)da. 
'OO 

In the energy estimates below, we make use of the fact that 

/oo -I pOO 

hd^fHd^f da^~- {[H,h]d^f)dc,f da s=: C||/j||h5/.+ 
OO -'J —oo 

Indeed, by integration by parts, 

/oo poo 
hdafHdaf da = - H{hdaf)daf da 
-OO J —OO 



/oc poo 
h{Hdaf)daf da - / ([i/, h]daf)daf da. 
OO J —oo 



Then, (fO]) follows by ([27]) and upon the observation that [H, h]daf = da{[H, h]f)- 
[H,dah]f. 
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It is readily seen that the energy funetion is equivalent to ||'u(<)||^3+3/2 + 

\\dtu{t)+u{t)dau{t)\\'jj, . Furthermore, ^'^{t) is equivalent to ||w(t)||^^+3/2 + ||9tu(t)|||^, 
if s > 1/2. 

We now state and prove the nonlinear energy estimate for (|6.ip . 

Proposition 6.1 (The nonlinear energy estimates). // (u, v) £ i7'*+'^/^(M) x H-''{R) 
for s > 1/2 solves (j6.ip on the interval t G [0, T] Jar some T > and if \\u\\^5/2+ < 
+00 for < t < T then 

(6.5) e'(t)<Ci forO<t<T 
and subsequently 

||u(t)||ff.+3/2 + \\dtu{t)\\H^ <C2 forO<t<T, 
where the constants C'i,C2 > depends on ||u(0)||^s+3/2 + ||c)fu(0)||H= • 
Remark 6.2 (Remark on the energy expression). One may try 

/•OO 

(6.6) / {{d^dtu)^ + d'^+\iHd^+\) da 



as the fc-th energy function for ()3.7|1 . Due to the multi-derivative nonlinear term 
u^&^u, however, the application to ()6.6p of the standard energy method is unwieldy. 
Indeed, one takes the i-derivative of the energy function and substitutes dfU by 
p.7p . but to arrive at an expression containing 



{didtu)di{u''dlu)da, 



which cannot be controlled by the energy function 

The idea of the proof is to write (|3.7p as a system (|6.ip to resolve the multi- 
derivative term u^&^u into two single-derivative terms udaU and udaV, which work 
favorably in the application of the energy method by canceling higher Sobolev 
norms when integrated by parts. 

Proof. We begin by investigating the time derivative of £^ by calculating 



(6.7) 



-£\t) = j J{dl+'dtu)Hd^{dl+'u) + [dldtv)[dlv)) da 



-■2 ■ 

Let us first compute E\. By using the first equation in ()6.ip and the integration 



by parts we may write 

/OO 
a^+^(-u9„w-i-w)i?aa(a^+iu) da 
-OO 

u(dl^^u){Hdl+^u) da 



(6.8) 



/OO 
(d^^^v)Hdo,d^^^u da + (lower order terms), 
-OO 



where (lower order terms) is a collection of terms which can be bounded in terms 
of energy in a routine way. The second inequality uses (|6.4p . Note that 

u{d^+''u){Hd^„+^u) da ^ \\u\\hW2^\\u\\1,+,. 
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Similarly, we compute 

/>oo 

^2 = / di{~udo.v + Hdlu + R{u, v)){d^v) da 

nOO 

(d^iud^v) - ud^^+^v){d^v) da ~ / {d^+^v)Hd^+^u da 



(6.9) 

+ / i9a''^)daRiu,v) da + (lower order terms). 

— oo 

Again, (lower order terms) is made up of terms which can be bounded in terms of 
the energy in a routine way. Note that 

\\d^{ud^v) - ud^+\,h2 ^ C\\u\\H4dav\\m-^- 

The third term on the right side of (j6.8p and the third term on the right side of 
(|6.9p cancel when added together. Other terms are bounded in terms of the energy, 
provided that ||u||//5/2+ < +oo. 

We have proved 

^ Cg^(l + £^)f, 
dt 

for some positive constant C and for some p > 1. The proof then is complete by 
applying Gronwall's inequality. □ 



Now, we make a few remarks on the existence, uniqueness and continuous de- 
pendence for p.7|) . or equivalently, for (|6.ip . In order to establish the existence 
of solutions, we would need to regularize the equation in a certain way. The most 
straightforward way is to introduce moUifiers (approximations to the Dirac delta 
function) into the right sides of (|6.ip . We then repeat the argument in the proof 
of Proposition 16.11 for the regularized problems to obtain energy estimates of the 
kind in (|6.5[) independent of the mollification parameter. Subsequently, the Picard 
theorem for ordinary differential equations on a Banach space applies to assert 
that solutions to the mollified equations exist on short intervals of time. The solu- 
tions can be continued on a time interval which is independent of the mollification 
parameter thanks to the uniform bound from the energy estimate. 

Next, an estimate similar to the energy estimate but in a low norm {{u,v) S 
i?^'(M) X i?^'-3/2(]R) with 2 ^ s' < 3 ) estabhshes that as the mollification param- 
eter tends to zero, the solutions of the mollified equations converge to a solution 
of the original (non-molified) equation (j6.ip . By interpolation, we find that this 
convergence occurs in the high Sobolev norms, as well. 

Uniqueness follows by the energy estimates for the difference of two solutions, 
and continuous dependence follows from the time-reversibility of the equations. 

The detail of local well-posedness via the argument above is carried out in [2] 
for the vortex sheet problem with surface tension. 

We summarize our result. 

Theorem 6.3 (The local well-posedness). The initial value problem for p.7p . pre- 
scribed with the initial conditions {uq,ui) £ H'^{R) x H-''^'^/^(R), s > 2 + 1/2, 
is well-posed on the time interval [0,T] for some T > and {u{t) , dtu(t)) e 
C([0,T];i7-^(R) X H''-^/^{R)). 
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Remark 6.4. The assumption that s > 2+1/2 is due to the commutator estimate 
used in the proof of the energy estimates ( Proposition \6.1] }. Although this commu- 
tator estimate can be improved by putting 1/2 derivatives on each of two copies of 
u appearing in the energy calculation, we need 2+ derivatives to couple with the 
transport equation. In this sense, the local well-posedness should hold at s > 2, 
although we do not prove it in this paper. 

7. StRICHARTZ estimates for the NONLINEAR PROBLEM 

At last, we are in a position to prove the Strichartz estimates (|1.3p and (II. 5|) for 
solutions to p.7p . 

Let us consider the initial value problem for (j3.7[) with the initial conditions 

u(0, a) = uo(a) and dtu{0, a) ~ ui{a), 

where (mo,mi) € H'{R) x i/''-3/2(M) for s > 2 + 1/2. Theorem [O] applies to 
ensure that a unique solution exists on the time interval [0, T] for some T > and 
{u{t),dtu{t)) G C{[0,T];H'{R) x H'-^/^{R)). 

We apply to (j3.7|) and the above initial conditions to obtain 

I (df ~ Hdlu + 2uda,dtu + u^d^d'^u = R{u, dtu), 
|i9„u(0, a) = 9^uo(a) and dtd^u{0,a) ^ d^ui{a), 

where 

(7.2) R{u,dtu) = diR+[dt,,u^dl + 2uda,dt]u. 

We view the initial value problem (|7.1[) as a linear problem for d^u of the form in 
(|5.ip and (|5.2p . where the solution u is the coefficient function V{t, a) and R{u, dtu) 
is the inhomogeneous term R{t, a). This can be done thanks to uniqueness of these 
initial value problems. 

Let us take s ^ 1 sufficiently large so that u e H\[0,T])H''{R), where > 1 
are large enough for the results in Section|4]and Section[5]and so that {d^UQ, d^ui) G 
L2(R) X H^^/'^{m.) and R G L'^{[0,T])Ll{R). Therefore, Theorem applies to 
yield the semiclassical Strichartz estimates (jl.Sp , and Corollary 15.31 apply to yield 
the fixed time Strichartz estimates (|1.3p . This completes the proof of Theorem ll.il 
and Theorem 11.21 

Appendix A. Local smoothing effect for the nonlinear problem 

Remarked here is how the method of positive commutator yields the local smooth- 
ing effect for p.7p of the gain of a 1/4 derivative's smoothness. The method used 
here was suggested to us by T. Alazard, N. Burq, and C. Zuily. 

Theorem A.l (The local smoothing effect). If s > 2 + 1/2 is sufficiently large, 
then the solution u G C([0, T]; II'^{R)) of the initial value problem for (|3.7p with the 
initial condition {uq,ui) G H^iR) x H^^'^/'^{R) satisfies the estimate 

(A.l) II {a)-" D:^+^/^uUHy,^T\)Li ^ C{\\u4h^, \\ui\\hs-v^), 

where p > 1/2. Here, (a) = (1 + a^)^/^ describes weighted Sobolev spaces and 

Da = -ida- 
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Proof. We recall that the energy method yields the local well-posedness of the initial 
value problem for dST]) and u e C([0, T];H^{R)) for some T > 0. 

Once we settle the issue of existence and uniqueness, we view (|3.7p as 

dfu - Hdlu + 2V{t, a)dtdau + V^{t, a)dlu = R{t, a), 

where the variable coefficient V{t^ a) replaces the solution u and R{t, a) replaces 
the remainder dtu), as in Section[7l We choose s > 1 large so that V{t, a) and 
i?(t, a) have the required regularity. 

Let ij} € C°°(M) be such that = 1 for |C| ^ 1 and = for |^| ^ 1/2, 
and let 

U = W + V, 

where w = ip(Da)u is the high frequency part of the solution and v is the low 
frequency part of the solution. It is straightforward that v belongs to every Sobolev 
space, and thus v is as smooth as we want. It suffices to show (|A.ip for w. 

Further, let ip{Da) = ii^iDa) + ip~{Da), where ip'^ = V'IfSso and ■0" = "^Ic^Oj 
and let = ip^{Da)u. We will present the argument only for positive high 
frequencies . For simplicity of notation, in what follows, we write w for w^. 

For positive high frequencies ^ > 0, the Hilbert transform H behaves like the 
multiplication by —i, and thus 

(A.2) Pw = {dl + idl + 2y(t, a)dtdc, + V^{t, a)dl)w = R{t, a), 

where 

R = ,p+{D^)R{t, a) + [2V{t, a)dtda + V^t, a)dl,^P+iD^)]u 
satisfies the estimate 

(A.3) ||i?||L=„(K) ^ U+{D^)R\\mR) + C{\\w\\hi^r) 

where C > depends on V and its derivatives. 
Let 



\dtw\\ 



A{a) = / (py^" df3 

J —OQ 

for p > 1/2 as in the statement of the theorem. Let ( 
L^-inner product. We compute 



•)^2 denote the (Hermitian) 



2i Im 



AR, w) dt^ 

L2 



(^AR,wj^^ dt - 



{APw,w)^2 dt 



w, AR) dt 
/l2 



{Aw,Pw)j^2 dt 



{[A,P]w,w)^2dt+ / {{P - P*)Aw,w)^2dt 



{Adtw, w) 



{Aw, dtw)^ 



Here, P* denotes the adjoint of P. The second equality uses that A is self-adjoint, 
and the third inequality uses integrations by parts in t. The sum of the boundary 
terms is 2i[lm {AdtW,w)^2]o- 

It is important to note that although the solution to the nonlinear equation (j3.7[) 
may be assumed real- valued, the solution if of (|A.2|) is not, since we have cut it off 
to positive frequencies. Indeed, if we add back in the negative frequencies so that 
w is once more real-valued and if we try to construct an appropriate commutant, 
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then we are inevitably led to use the Hilbert transform as part of the commutant, 
and the resulting boundary terms do not cancel. 
Next, it is straightforward that 

[P, A]w ~ SiAad^W + SiAaadaW + iAaaaW 

+ 2V{t, a)Aadtw + 2V^{t, a)AadaW + V^{t, a)AaaW. 

r r 

Now, we compute / {Aw, Pw) dt = / (P*y4.w, ?i')^2 fit. Integrations by parts 
Jo Jo 

in t and in a yield 



{Aw,d1w) ^.^dt ^ {Aw , dtw) — {Adtw,w)j^2 + f {Ad^w,w) ^^dt 



2i 



lm{Adtw,w)^2 + / (^Ad^w,w) dt 



T j-T 

(^Aw,id^w) dt = / {ida{Aw),w) j^2 dt. 
Jo 



Similarly, 



{Aw, 2Vdtdaw)^2 dt = {2VAw, 9„ 



/L2 



{2VtAw,daw)^2 dt 



+ j (^{2Va.Adtw,w)^2dt+{2Vdtdc.{Aw),w)^2)dt, 



{Aw,V^dlw)^,dt= [ (^{{V^)^^Aw,w)^, 



+ {2{V^)MAw),w)^2 + {V^dl{Aw),w)^2)dt. 



Therefore, 



1 T 





3i / {AadaW , daw) j^2 dt ~2i J lm(^AR,wJ dt + Ai lra{Adtw,w) 

+ i J (^{AaaaW,w)]^2 + {2VAadtW,w)]^2^dt 

+ j (^{2V'^Aa,daW,w)^^ + {V'^A^a,w,w)^.,^dt 



{2VAw, daw) 



{2VtAw,daw)^2 dt 



+ J (^{2Vo,dtiAw),w)^2 

+ {{V^) o^a Aw, w) ^2 + {2{V^)MAw),w)^2)dt. 
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Since A and its derivatives are in L°°(M), using the energy estimates, Sobolev 
embeddings, and the estimates for R, it follows that 

Ucd^Mlh^Li ^ ^^(Plli|H-3/4 + lklli|H3/4 + Ml^HV^ + \\dtw\\l^H-v^ 

+ (1 + \\VYl-hO\M%^h^J2 + (1 + \\V\\l^H,)\\w\\l^^.j, 
+ {l + \\dtV\\l^H.)\\w\\l,^^,j, 

+ (l + l|a^(^')lli?.LS=)ll^lli|L2 

^ C(r, ni^3/.^,)(lkll,o,^3/4 + ||a,Z.||^^^-3/4f 

for some s' > 0. This completes the proof. □ 

Appendix B. Assorted proofs: formulation 

We collect the proofs of p.2ip . Corollarv l2.5l Lemma [321 £^nd Lemma l373l 
In order to estimate ||7t||Hs in terms of u and 9, we recall that 

da = sdle + a„((c/ii - w • t)7) - 2Wt • t - ^79^7 + 2(c/ii - w • t)w„ ■ t. 

We expand • t as 
Wt-t=Re(z„$(WO) 

-^"(2^-"(")P^/ .M-"?K) '"') 

+ Re(^..(.)Pv/7(a')(5^ff^^a') '-^ JM,. + 

where J^[z]f ~ Re (^2iZa,H (^j-^ + Zaia)lC[z]f^ . Accordingly, the above equation 
for 7( takes the form as 

(B.l) (l+2^[z])7t = 5a20+a„((C/"-W.t)7)-i7a„7+2(C/ll-W.t)W„.t-2i?5. 

It is proved in [2, Lemma 6.1] that (1 + 2J'[z])~^ : ^ is bounded. One 
observes that i?5 may be written in such a way that it is the sum of terms which 
differ from R3 and R4 by multiplication by i, and therefore, they are estimated 
mutandis mutandi to yield that 

for s ^ 1. The argument in the proof of [2, Lemma A. 4] applies to assert (|2.21[) 
Proof of Corollary \2.5l We write 

Here, the last term is related to R3 and R4 in the proof of Lemma [^T^ and thus it 
is estimated in a similar way. 
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By the usual product rule, (|2.24|) follows. 

Proof of Lemma W^ For s = 0, 1, we use the transport equation 
plication by 9 to (|3.8p and integration by parts yield 



□ 

By multi- 



^ / e^da = i / e'^dau da 



j 9Hdauda + J 9ri{t,a) da. 



Wc obtain the analogous identity for J and by adding, 
d 



dr 



Gronwall's inequality then applies to give that 



0{t)\\H^ ^II^oIIhi + / {\\u\\H- + \\r,\\H^)dt' 



C 



\do.u\\L^ (^\94h^ + ^* {\\u\\h2 + llrillffoj exp (^j^^ Wd^W 
^C{\\u\\H2)(l+\\u\\H^ + \\r,\\Hi). 

Indeed, Wr^u^ Cm\H-){^ + Mm). 

Next, for s = 2 by multiplying (|2.10b ) by d^9 and by integrating it it follows 
that 

\\dle\\l. ^ WdtuU. + \\u\\\, + wwl^WdML^ + cmm){^ + Mmf. 

Together with the 116*11^1 estimate above, this proves p.9p for s ~Q. For s > 2, we 
take derivative of (|2.10b ) and repeat the argument. This proves the assertion. □ 

Proof of Lemma lKR First, it is straightforward to see that 

dtri = ~H{mt ■ t) - H{m ■ n)et + mt • n + (m • t)et, 

where 



~2i 



H, — 



7a 



+ Za_dt ( /C[z 



2z^* 



la 



Then, (P?^ and (P?^ apply to assert that 

WmtWrn ^ Cmu\\mA\u\\H^^^, \\^tu\\Hs-^), 

and, in turn, it follows p.lOp . The difference is estimated in the usual way. 
Next is the estimate for dtr2. We recall from the proof of Lemma ^TM that 

dtr2 = \Hdh + dt{Ri + R2 + R3 + Ri). 
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In order to estimate for 9^ 7, we take the derivative with respect to t-variable of 
(IB. II) to obtain 



(td + J[z])jtt = - Re( 2^z„ti^ ( — ) - 2iz^H ( ^z^t 



r-\ ^ Zarj fit \ Za f , zt{a) - Z((a') a 
+ 0aat + 5ta„(7(C/ll -W.t))-29t fi77„-(;7ll -W.t)W„.t 



2< 



^4 

Each term on the right side of the equation is estimated by using various estimates 
we estabhshed previously, and then we conclude that 

\\dh\\H^ ^Ci\\u\\Hs+2,\\dtu\\H. + i). 

Again, using the estimates established previously, we obtain 

\\dt{Ri + R2 + R3 + R4)\\h^ ^ Ci\\dtu\\Hi,\H\H^+^A\dtu\\H^-i). 

The differences of dfj and dtRj's, j — 1,2,3,4, are obtained in the usual way. 
Therefore follows ^J^. 

That is, without the cancellation of the highest-order derivative term daOd^u in 
r2, the remainder R(u,dtu) is of second-order in a. This completes the proof. □ 

Appendix C. Assorted proofs: parametrix construction 
We first derive (|4.10p . Let us first write iHD^w via the Fourier transform as 

iHDlw{t,a) = -Hdlw{t,a) 



J_ // e^(a-"')e|^|3 
47r 



In what follows, we recall that wc implicitly assume that both ^ and ipa are large 
and positive and that ^ ~ 2^ . Our goal is to eliminate the dependence on ^ in the 
above integral so that the integration in a' and ^ yields 6 {a — a') (so that the above 
representation reduces to an integral in f3 and ^' only). To this end, we write 

tp{t, a, C') = fit, a',C) + ^(a, - 
and we perform a change of variables to obtain 

iHDlw{t,a) = ^ jj e^("-"')''|77-H^>(a,a')|^ 

We further write a') = (pa{t, S,') + Oi, a' , ^'), where 
<Pi{t, a, a', = ]^faa{t, a', - 
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for some <1> = 0(supjj I^Qipl), and accordingly, the above integral becomes 
iHDlw{t,a) 



47r 



We keep in mind that ipa in the above expression is evaluated at {t,a',£,'). Now, 
'Pi is a smTL of terms multiplied with powers of a — a', which upon integrations by 
parts in rj are cancelled and the above integral, in turn, becomes 

(C.l) '^^^ \„ ^ 

• // e-'^^' e''^^'-°''^'\f{p) df3dC da'dr,. 

Under the assumption that either both rj and tpa arc large and positive or both 
are large and negative the above formal argument is justified. Indeed, the dyadic 
frequency localization assumption (|4.2|) on / implies that is also localized to 
dyadic frequencies (possibly with different constants), and hence the singularity of 
|?7| at ?7 ~ docs not enter into the above calculation. 

We now expand [77 + ipa\^ for both 77 and large and positive (see Lemma 14.81 
for a justification of this) 

\ri + ^al' = hi' + 3|^„||77|2 + 3|^„nr;| + |^ap. 

Again, we keep in mind that ipa is evaluated at (t, a' ,^'). Substituting this in (jC.ip 
and integrations by parts in a' then yield that 

iHDlw{t,a) = ^ JJ e'^"^"')" 

• (^|77p + 3|77p|(^a| + 3|77||(^Qp + \ipa\^ + 3i{\r]\ + {ipaD^Paa " ^aac^j 

;-*^«'e"^(*'"'«')/(/3)d/3dC'da'd77. 



Finally, (|4.10p follows upon integrations in a' and rj. 

Next, we show that the mapping in (|4.20[) is invertible for each C G [— e, e] and 
5$ t ^ 2^/'^T. It suffices to show that 



da 



^ C"^ > for t 



8(3 

First, if rjit) is a solution of (j4.19b ) with the initial condition 77(0) = C then 
[(1 + 2-^'^e'^vf\ = 2(1 + 2-^'^e'''v)'i'''H^'''il = 2K.(t, 2^'^a){l + 2-^/2^^/2^)2^ 
whence 

-c(i + 2-^/^e''^vf ^ [(1 + 2-='"~e'^if\ ^ c{i + 2-=/-^^'^^? 

for some C ^ ||V'Q|jL^Lc>= . We recall that the dot represents differentiation in the 
f-variable. By Gronwall's inequality it then follows that 

(1 + 2-^/2^1/2^)2 cxp(_Ci) ^ (1 + 2-''^e'^vf < (1 + 2-'^'^e'^Cf cxp(Ci). 
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That is, (1 + 2-J'/2^i/2^(i))2 ^ _^ 2-J'/2,e^/2()2(l + Oit)). 
Next, we calculate 
9(1 + 2--''/2^i/2^) 



dp 

^2 0-7/2.1/2 ^(1 + 2--''/^^^/^^) dr, 
^ dp dp 



' dp 



d{l + 2-J/2gi/2?7) 
dp 



2^/V„„(t,2^"/2^)(l + 2-^/Y/2r;) 



da 
dp 



< (2|14| + 2'^'\VU1 + 2--'"/Y/^)|') 

By Gronwall's inequality it follows that 
9(1 + 2-^/2^1/27;) 



9(1 + 2-^/2^1/2^) 
dp 



dp 

for some C, C" > 0. 
Finally, we calculate 



(i)<exp(Ci2^/2||V^||2^^.,^)(C' + 2^/ 



* fda 
\dP 



r/Sax 2 



2 9q; 96; 
9/39^ 



9a/3 ^^./2 1 
5/3 4' ^ 



/2(1 + 2-^V2^1/2^)-l/2 5(l + 2^^^?^ _ 2^"/2^-l/2T4(t, 2^/2c 

9p 



9c 



> _(2.V2^-i/2|y^| + 1) ( ^ ) - C(l + 2-^/2^1/2,?) 



^ 9(1 + 2-^/2^1/2^) 
9/3 



9c 



9/3 



-C sup 2-^/2 



9c 



dp 



The claim then follows by Gronwall's inequality once ^ 

2~i/'^T are observed. Consequently, the inverse function theorem applies to give 
that the mapping ([4201) is invcrtiblc for < i < 2-^/'^T. 



1 and ^ t ^ 



Appendix D. Energy estimates for the linearized equation 

This appendix concerns the energy estimates (|5.3[) and (|5.4[) for the linear prob- 
lems (|5.ip and ()5.2|) , respectively, the idea of which will be used repeatedly through- 
out this work. 

We write the second-order equation ()5.ip as the following first-order system 
(D 1) \dtu^ -V{t,a)daU + v, 

^ [dtv = ~-V{t,a)dc,v + Hdlu + Vtit,a)do,u + V{t,a)Vo,{t,a)daU. 

In other words, v ~ dfU + V{t, a)daU is the directional derivative. 
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Let US define the fc-th energy associated to the above system by 
(D.2) S'^it) = - / {{d',+'u)Hd.{d^+'u) + (9^)2) da 



2 . _ 

and define the energy function of order s by 

(D.3) e^^t) = \\u\\Ut) + \\v\\Ut) + J2£''it). 

fc=l 

Note that Hda is a positive operator witli the symbol of its Fourier transform |^| 
and that 

II 



= / {.f + IHd^f) da. 



It is immediate that ^^{t) is equivalent to ||^,_|_3/2 + ||9ttt(t) provided that 
bounded. 

Dy investigating the time derivative of by calculating 

' POO 

-£\t) = / {{di+^dtu)Hdo.[dt'^^u) + {d%v){dlv)) da 



(D.4) 



dt 



— OO 



The first equality uses that Hda is self- adjoint. 

Let us first compute . By using the first equation in (jP.ip and the integration 
by parts we may write 



5^= / di+^{~V{t,a)daU + v)Hda{d'^+^u) da 
Vit,a){d^+'^u)Hd^+'^u da 

/OO 
{d^^'^v)Hdad^^'^u da + (lower order terms) 
-OO 

/OO 
{d^^^v)Hdo,d^^u da + (lower order terms), 
-OO 

where (lower order terms) is a collection of terms which can be bounded in terms 
of energy in a routine way. The third equality uses (j6.4p . 
Similarly, we compute 

/OO 
di{~V{t,a)daV + Hdlu 
-OO 

- Vt{t, a)daU - y(t, a)Va{t, a)dau){d^v) da 

/OO 
{d^{V{t,a)dav) - V{t,a)dt'v){d^v) da 
-OO 

/OO 
{d^^v)Hd^'^u da + (lower order terms). 
-OO 

Again, (lower order terms) is made up of terms which can be bounded in terms 
of the energy in a routine way. The first term on the right side is bounded by 

\\V\\H4^c.v\\H>^-^. 
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The first term on the right side of 8^ and the second term on the right side of 
£2 cancel when added together. Other terms are bounded in terms of the energy. 
Therefore we have proved 

for some positive constant C, provided that ||t^||Loo([o.T])iy= °=(K) is bounded. The 
energy estimates (|5.3p with s ^ then follows once Gronwall's inequality applies. 
For the inhomogeneous problem (j5.2p with s ^ 0, we proceed similarly to obtain 

— <Ce:^ + ||i?(i)||?f.(K), 

from which (|5.4p follows. 

We next show that this can be extended to negative indices s. Let ^ s' ^ 3/2, 

and observe that if u solves (jS.ip for t S [0, Ti], then U = (ZJq)"* u satisfies 

{PU =[P,{D^)-'']u for <G[0,ri], 

\c/(0,a) = {Do)''' dtUiO,a) = {DoV' m. 

Thus U satisfies the energy estimates (|5.3p - (|5.4p for s ^ 0: 

(D.5) WdtUU 

^CTAm{Do.r]u\\L 2([0,Tl])ifJ + ||uo||^a-.'+3/2 + lluillj:^,-,')- 

The commutator term is bounded by 



s' + l, 



\\\P, (Da) ' ]w||l2([0,Ti])HJ C'(ll^t"llL2([0,Ti])Hr = ' + """^^([O.Ti])//^ 

^ c(ri)^/^(ll9*«IU.([o,T,])«r»' + ll"llL==([o,T,])«r='+^)' 

but for Ti > sufficiently small and fixed, this can be absorbed into the left hand 

side of p.5|) to get 

(D.6) 

l|5tM||^„c([o^Ti])ffr = ' + ll"IL~([0,Ti])ffr°'+'^' ^ Cdluoll^j-.'+a/s + )■ 

Applying this argument over a finite number of time steps of size Ti > yields the 
estimate for any finite T > 0, with constants dependent on T. A similar argument 
applies in the case of (|5.2p to yield the estimate (|5.4p for —3/2 ^ s ^ 0. 
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